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Time-dependent expectation values and correlation functions for many-body quantum systems are 
evaluated by means of a unified variational principle. It optimizes a generating functional depending 
on sources associated with the observables of interest. It is built by imposing through Lagrange 
multipliers constraints that account for the initial state (at equilibrium or off equilibrium) and for 
the backward Heisenberg evolution of the observables. The trial objects are respectively akin to 
a density operator and to an operator involving the observables of interest and the sources. We 
work out here the case where trial spaces constitute Lie groups. This choice reduces the original 
degrees of freedom to those of the underlying Lie algebra, consisting of simple observables; the 
resulting objects are labeled by the indices of a basis of this algebra. Explicit results are obtained 
by expanding in powers of the sources. Zeroth and first orders provide thermodynamic quantities 
and expectation values in the form of mean-field approximations, with dynamical equations having a 
classical Lie-Poisson structure. At second order, the variational expression for two-time correlation 
functions separates—as does its exact counterpart—the approximate dynamics of the observables 
from the approximate correlations in the initial state. Two building blocks are involved: (i) a 
commutation matrix which stems from the structure constants of the Lie algebra; and (ii) the 
second-derivative matrix of a free-energy function. The diagonalization of both matrices, required 
for practical calculations, is worked out, in a way analogous to the standard RPA. The ensuing 
structure of the variational formulae is the same as for a system of non-interacting bosons (or of 
harmonic oscillators) plus, at non-zero temperature, classical gaussian variables. This property 
is explained by mapping the original Lie algebra onto a simpler Lie algebra. The results, valid 
for any trial Lie group, fulfill consistency properties and encompass several special cases: linear 
responses, static and time-dependent fluctuations, zero- and high-temperature limits, static and 
dynamic stability of small deviations. 
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I. INTRODUCTION 

Variational methods have proved their flexibility and efficiency in many domains of physics, chemistry and applied 
mathematics, in particular when no small parameter allows perturbative approaches. In physics one may wish, for 
systems of fermions, bosons or spins, to evaluate variationally various quantities, such as thermodynamic properties, 
expectation values, fluctuations or correlation functions of some observables of interest. One may deal with the ground 
state, with equilibrium at finite temperature, or with time-dependences in off-equilibrium situations. 

We want, moreover, to perform these evaluations consistently. Suppose for instance that we have optimized the 
free energy of a system by determining variationally its approximate state within some trial class; nothing tells us 
that this state is also suited to a consistent evaluation of other properties. Is it possible, remaining in the same trial 
class, to optimize some other quantity than the free energy, for instance a statistical fluctuation? 

The wanted properties can be of different types. For instance, one may be interested in both the expectation 
values and the correlations of some set of basic observables. Consistency then requires the simultaneous optimization 
of these quantities. In order to evaluate them in a unified framework it appears natural, as in probability theory 
and statistical mechanics, to rely on characteristic functions, or more generally for time-dependent problems, on 
functionals that generate correlation functions. As in field theory, time-dependent sources ^j(t) are associated with 
the basic observables Qj. Expansion of the generating functional in powers of these sources will supply expectation 
values and correlation functions of the set Qj in a consistent fashion. 

Our strategy, therefore, will be the variational optimization of a generating functional V’lC} foi' (connected) corre¬ 
lation functions. To face this problem we will rely on a general method allowing the systematic construction of 
a variational principle that optimizes some wanted quantity. In this procedure, the equations that characterize this 
quantity are implemented as constraints by means of Lagrange multipliers. 

The desired generating functional is expressed as V’{C} = InTrA(ti) D in terms of two entities, the state D in the 

/ OO _ 

dt'''^f,j{t')Qf{t',t)] taken at the initial time 
i 

t = ti. The operators Qj{t',t) entering A{t) are the observables of interest in the Heisenberg picture, t' being the 
running time and t the reference time at which (5^(t,<) reduces to the observable Qj in the Schrodinger picture. The 
variational determination of together with its expansion at successive orders in the sources f,j{t') provides the 
various desired outcomes: namely, at zeroth order (for A(t) = I) thermodynamic potentials if D is a (non-normalized) 
Gibbs state, or ground state energy in the zero-temperature limit; at first order expectation values; at second order 
correlation functions for an initial off-equilibrium state D, etc... Static correlations within the state D will be obtained 
from sources located at the origin ti of times. 

Implementing the variational principle requires the use of formally simple equations which characterize the two 
ingredients D and A(ti) of the generating functional '0{^}, and which will be taken as constraints accounted for by 
Lagrange multipliers. The state D = exp(—/lAT) will be characterized by Bloch’s equation for D{t) = exp(—rAT); the 
associated Lagrange multiplier is an operator depending on an imaginary time r. In order to characterize the second 
ingredient A{ti) of V'{C}) have defined A{t) for an arbitrary initial time t. The observables Qj{t',t) entering the 
’’generating operator” A{t) then satisfy a backward Heisenberg eguation [Eq. (HOI) ] in terms of the reference time t 
which will eventually be fixed at A. This backward Heisenberg equation plays a crucial role as it produces for A{t) 
the formally simple differential equation CD; the associated Lagrange multiplier is a time-dependent operator. The 
equations for the density operator D{t) and for the generating operator A{f) are complemented by the boundary 
conditions D(0) = I and A(-|-oo) = I, where I is the unit operator. The ’’time” t of D{t) varies forward, the time t 
of A{t) backward. 

Unrestricted variations of the trial operators T>{t) and A{t), and of their associated multipliers, yield the exact 
generating functional, the stationarity conditions being the exact dynamical equations for D(t) and A{t). These 
operators should be restricted within a trial subspace to make the evaluations feasible, and then the resulting equations 
become coupled. Their solution will be simplified by expansion in powers of the sources f,j(t), yielding a tractable, 
unified and consistent treatment. 

In this article, we choose as trial subspace a Lie group of operators. This will be the sole approximation. The 
formalism will be developed for an arbitrary Lie group (^Secs. lHlllVl) . Explicit calculations are then allowed for a 
sufficiently simple underlying Lie algebra. 

For arbitrary systems and for any trial Lie group, mean-field like approximations are recovered at zeroth and first 
orders in the sources fSec. lIVp for thermodynamic quantities and static or dynamic expectation values (for instance, 
selecting for fermions at finite temperature the Lie algebra of single-particle operators yields back in this general 
framework the static and time-dependent Hartree-Fock theories). 

However, at second order in the sources, the formalism generates non trivial results for fluctuations, for static 
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correlations and for two-time correlation functions (Sec.|Vl). Remarkably, new variational approximations come out 
for these quantities, although the trial operators belong to the same simple class as the one that provides standard 
results for the expectation values (exponentials of single-particle operators in the example of fermions). Within the 
Lie group the trial operators adapt themselves to each question being asked so as to optimize the answer - while the 
use of the generating functional ensures the consistency of the results thus obtained. 

The second half of this article ISecs. lVUlXl is devoted to the properties of the outcomes of the present variational 
theory with trial Lie groups. In particular, the static and dynamic stabilities are related to each other (Sec. IXl); a 
quasi-bosonic structure fSec. IVlUl) and a classical structure (Sec.|2) are exhibited for any system; various consistency 
properties are reviewed (Sec. ED). 

Some past works are related to the present variational approach, which generalizes and unifies them within a natural 
framework. For fermions at finite temperature, the optimization of expectation values has been shown to lead to the 
static HF and dynamic TDHF equations Q , while variational expressions for fluctuations Q and correlation functions 
Q have been derived. In p articular, the large fluctuations observed in heavy ion nuclear reactions have thus been 
correctly reproduced p-RlLas recalled in Sec. IVI ^plications of the variational principle of Sec. lIICI have been 
made to bosonic systems |l2l|. to Bose condensation [ij, HI, to field theory in including two-time correlation 
functions [l^ . [l^ , and to restoration of broken particle-number invariance for paired fermions at finite temperature 
Let us also mention an application to cosmology [l^ and an extension to control theory [I^ . In the case 

of small fluctuations around the mean-field trajectory for fermions, the variational principle leads to time-dependent 
RPA corrections similarly to other approaches assuming either initial sampling of quantum fluctuations [2ll . l22j| or 
directly solving time-dependent RPA (for a recent review, see M)- 

The main results are recapitulated in Sec. m before the conclusion in Sec. IXIIl 


II. A UNIFIED VARIATIONAL APPROACH 
A. Generating functional 

We consider a quantum physical system prepared at the initial time ti in a state represented by a density operator 
of the form D oc exp(—/SAT) in the Hilbert space The dynamics is generated by the Hamiltonian H, possibly 
time-dependent. If the system is in canonical (or grand canonical) equilibrium, one has K = H (or K = H — fiN); 
for dynamical problems K needs not commute with H. Ground state problems are trailed by letting /3 —>■ oo. 
Partition functions will be evaluated as Tr D for the unnormalized state 

D = e-^^. ( 1 ) 

The normalized state will be denoted as D = D/TyD. We are mainly interested in expectation values, fluctuations 
and correlation functions of some set of observables denoted as Q® in the Schrodinger picture. We will work in the 
Heisenberg picture in which the observables 

Qf{ti,ti) = U^{tuU)Q^U{tf,tO ( 2 ) 

depend on two times, the initial reference time ti and the final running time t{. In the unitary evolution operator 

U{tf, ti) = Te~^ (3) 


T denotes time ordering from right to left. 

In order to generate consistently the desired quantities, we associate with each observable Qj a time-dependent 
source (t) and we introduce the generating operator 




( 4 ) 


Then, the generating functional 


V>U} = lnTrA(fi)D, (5) 

which depends on the functions fj{t), encompasses the quantities of interest. In particular the partition function 
Tr D, the expectation values 


{Qj)t = TTQf{t,ti)D 


( 6 ) 
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at the time t, and the two-time correlation functions 

Cjk{t',t")=TrTQfit',U)Qf{t",U)D-{Q,)t'{Qk)t'', (7) 

are obtained as functional derivatives with respect to the sources according to the successive terms of the expansion 
of 

V'{^} = In Tr D -P f /“ dt 

-I ItT dt' r dt" ik{t") CMt',t") + ■■■■ 

Variances AQ‘j{t), hence fluctuations AQj{t), are found as Cjj{t,t)-, static expectation values and correlations in the 
state D are found by letting both times equal to tj. Linear responses are also covered by the formalism. 


B. The constraints 


In order to optimize simultaneously all the quantities embedded in the generating functional use is made of 

a general procedure Qi inspired by the Lagrange multiplier method. The purpose is to construct an expression 
whose stationary value provides the quantity we are looking for, namely here V'{^} = InTr A(ti)ZI. To implement in 
this expression the quantities D = and A{t), we will characterize them by equations regarded as constraints. 

To characterize the state D, we introduce a trial ’’time”-dependent operator 27 ( t ) compelled to satisfy the initial 
condition 27(0) = / and the Bloch equation 


^^ + KV{t)=0, (9) 

dr 

where the imaginary ’’time” r runs from 0 to /3. One recovers D = exp(—/I2L) from the solution of ([S]) for t = (3. 

Let us turn to the generating operator A{t) defined by (jd]) in terms of the observables In this form, A{t) 

looks difficult to handle and we wish to characterize it by a formally simple equation that can be taken as a constraint. 
To this aim, the operators t) are regarded as functions of the initial time t rather than of the final running time 

t' which is kept fixed at t' = tf. They thus satisfy the backward Heisenberg equation Q 


dQ 




( 10 ) 


This differential equation, together with its final boundary condition Q^{tiAf) = Q^{ti) at t = tf, is equivalent to 
the definition ([2]) of Qj{tf,t). Contrary to the standard forward Heisenberg equation (a differential equation in 
terms of t'), the backward equation (fTUl) holds even when the observable Q® and/or the Hamiltonian H depend on 
time in the Schrodinger picture. Note that in the backward Eq. m, H is written in the Schrodinger picture if it is 
time-dependent. 

In the present context, the forward equation for Qj{t',t) would be of no help in dealing with the definition ([4]) of 
the generating operator A{t) whereas the backward Heisenberg equation (ITUl) readily provides the differential equation 


^ [A{t), H] + tAit) ^ c,(t) gs = 0, 

3 


( 11 ) 


which, together with the boundary condition H(-|-oo) = J, is equivalent to (|4]). The generating functional'!/’{'?} defined 
by ([5]) involves Aitf), and this operator will be found by letting t run backward in Eg. dill) from -|-oo to with the 
final condition H(-|-oo) = I. Here again we shall simulate the operator A{t), solution of the exact equation (fTT|l . by a 
trial operator A{t) which will satisfy dill) approximately. 

In order to optimize the generating functional '0{^} = InTrH(ti) D, a variational expression [Eq. (|T^ below] will 
be constructed, which relies on the equations dH]) for D and (fTTl) for A{t). These equations are regarded as constraints 
with which Lagrange multipliers will be associated. We denote the Lagrange multiplier accounting for the Bloch 
equation (jH]) by A(t), an operator depending on the ’’time” t ; we denote the Lagrange multiplier accounting for 
Ea. dTTI) by V{t), a time-dependent operator. These notations are inspired by the duality between observables A and 
states 27 at the root of the algebraic formulation of quantum mechanics where expectation values are expressed 

as scalar products Tr^27. Here, 27(r) (for 0 < r < /3) and the multiplier 27(t) (for t > t,) are state-like quantities, 
whereas the multiplier A{t) and the operator A{t) are observable-like quantities. 





6 


C. The variational principle 


The implementation of the constraint ([9]) for 'D(t) by introducing the Lagrange multiplier A{t), and of the constraint 
(fTTI) for A{t) by introducing of the Lagrange multiplier T’(t), results in the variational expression @ 


'i>{A,V} = InTrAit = ti)V{T = / 3 ) 
'dV{T) 


( 12 ) 


dr Tryl(T) 


dt Tr 


dr 

dA{t) 1 
dt if 


+ KVir) 


[Tr A(t)V{t)]- 


[A{t), H]+iA{t)Y.^,{t)Q' 


V{t)[Tv A{t)V{t)] 


-1 


where normalizing denominators are included for convenience. Together with the mixed boundary conditions 


A{t = +oo) = I, Vij = 0) = / , (13) 

'^{A.V} should be made stationary with respect to the four time-dependent operators T’(t),^( r),^(t),T>(t) (with 
0 < T < l3 and U < t < -|-c»). The stationarity conditions include the additional continuity relations 


'D{t = fj) = 'D{t = ti) and A{t = (3) = A{t = t {), 


( 14 ) 


another argument for the notation. (In view of this continuity one might replace t by a complex time t = ti + i{(3 — T)h 
so as to rewrite the two integrals of (HU) as a single integral on a Keldysh-like contour @,[ 23 .) 

For unrestricted variations of A and T) the stationary value of dr is the required generating functional '0{^}. It 
is attained for values of A and V that let the two square brackets of m vanish, so that we recover the evolution 
equations and CD, the solutions of which are T’(t) = exp {—tK) and A{t) = A{t). 

The data of the problem, K,H and the observables Qfj are operators in the Hilbert space Jif. They all enter 
explicitly the variational expression ’^{A,'D'\. Simpler variational principles (VPs) derive from (1121) in two special 
circumstances. If the initial state D is workable, the first integral over r should be omitted. In this case, the variational 
principle (fT^ can be viewed for ^j{t) = 0 as a transposition of the Lippmann-Schwinger VP [l^l from the Hilbert 
space, with duality between bras and kets, to the Liouville space, with duality between observables and states [1]. 
For static problems, the last integral over t is irrelevant [13. Classical problems enter the same framework, with the 
replacement of the Hilbert space by the phase space, traces by integrals and commutators by Poisson brackets. 

As usual, practical exploitation of the above variational approach relies on restricting the trial spaces so that the 
expression CD of '^{A,'D} can be explicitly worked out. The denominators Tr AV have been introduced so as to let 
the functional 'T{A^'D} be invariant under time-dependent changes of normalization of T) and A. This allows us to 
select a ’’gauge”, that is, to fix at each time the traces of T) and A in such a way that the stationarity conditions take 
the form 


Tr SA{t) 
dA{T) 


Tr 


dr 


dV{T) 
dr 

— A{t)K 


KV{t) 


= 0 , 


5V(t) = 0, 


(0 < r < /3) 

(0 < T < /3) 


Tr 


^ 1 [A{t), H] + ^Ait) 

Tr 6A(t) 




dt ih 


5V{t) = 0, (t > U) 

= 0, (t > ti) 


(15) 

(16) 

(17) 

(18) 


in the restricted space for A{t)^ A{t), T>{t) and the corresponding space for their infinitesimal variations. In 

agreement with the boundary conditions (|T^ and (|TD relating the sectors r and t, equations IT^ and (|T8l) for V 
should be solved forward in time, with r running from 0 to /3 and t running from t, to 00 , whereas Eos. (1171) and (1161) 
for A should be solved backward. We obtain the stationary value of dt as 


tf{^} ~ In Tryl(t) 'D{t) = In Tr^(T) 27 ( t ) ( 19 ) 

for arbitrary t or r, as a consequence of the stationarity conditions written for 6A oc A, ST) oc 27. In the following the 
restricted choice of the trial space will imply that the allowed variations 5A around A depend on A, and likewise for 
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(527, so that the forward or backward equations are coupled. Practical solutions will take advantage of their 

expansion in powers of the sources 

The variational procedure has duplicated the dynamical equations, introducing Eqs. (HH) and (113, besides the 
approximate Bloch equation m and the approximate equation (fT71) for the generating operator A{t). While the 
formalism was set up in the Heisenberg picture, the stationarity condition (1181) reduces, in the absence of sources and 
for unrestricted variations of A(t) and 27(t), to the Liouville-von Neumann equation 


dV{t) 

dt 




( 20 ) 


The Lagrange multiplier matrix T>{t) thus behaves as a time-dependent density operator in the Schrodinger picture. 
However, such an interpretation does not hold in the presence of sources, in which case 'D(f) is not even hermitian. 


III. LIE GROUP AS TRIAL SPACE 
A. Parametrizations and entropy 


From now on, we specialize the trial space for the operators A and 27 involved in the variational principle (I12L 
assuming it to be endowed with a Lie group structure. This will be our sole approximation. The Lie group is generated 
by a Lie algebra {M} of operators acting on the Hilbert space a basis of which is denoted as M^. This algebra is 
characterized by the structure constants entering the commutation relations 

( 21 ) 

(We use throughout the convention of summation over repeated indices.) These constants are antisymmetric and 
satisfy the Jacobi identity 

It is convenient to include in the algebra the unit operator I, denoted as Mp. A seminal example is provided, for a 
many-fermion problem, by the algebra of fermionic single-particle operators = cifcii, [with a = {v^p) 0]. Other 

examples are given in Sec. lXI A1 such as, for condensed bosons, the set of creation and annihilation operators and of 
their products in pairs. 

The operators A and 27 are then parametrized, at each time r or t, according to 

(-22) 

The parameters L°‘ and J“ are functions of the times r or t; their sets will be denoted as {L} and {J}. They will 
be complex due to the presence of the sources £,j{t), and to a possible non-hermiticity of the operators M^. [In 
the fermionic example, the operators 27 have the nature of non-normalized independent-particle states; for bosonic 
problems they would encompass coherent states.] 

The results will be conveniently expressed by writing 27 as 27 = Z27, where Z denotes the normalization factor 

Z{J} = Tr27 = Tre“'“'^“ , (23) 


and where 27 is a normalized operator. Instead of the set {J}, the operator 27 can alternatively be parametrized by 
Z and by the set {2?}, defined by 

iJlu Z 

2?„=TrM„P=^ (24) 

and including Rq = 1. [In the example of fermions, for a = {v^p) 7 ^ 0, the set Ra = Ru^i = Trala,y27 are the Wick 
contractions associated with the independent-particle trial operator 27, so that a covariant vector {2?} can also be 
regarded as a single-particle density matrix. Contravariant vectors such as {L} or {J} are then regarded as matrices 
with switched indices, so as to produce the usual expressions in Hilbert space for operators and scalars, such as 
L“M« = or J“2?a = 

The converse equations of (1241) . which relate the set {J} to the set {2?},Z, involve the von Neumann entropy 

{kn = 1 ) 


Tr 27 In 27 


= InZ- J“ Ra , 


S'{2?} = -Tr 27 In 27 = lnTr27 


Tr27 


(25) 







(26) 


J“ = 


dRa 


(a ^ 0), 


J° 


\nZ - S{R}-^r Ra- 

a/O 


The Legendre transform (I^5l) - (I^S1) from lnZ{J} to the von Neumann entropy 5'{i?} stems from the exponential 
form of T) as function of the parameters {J}. The situation is the same as in thermodynamics where a Legendre 
transform relates thermodynamic potentials and entropy when going from intensive to extensive variables, due to the 
the Boltzmann-Gibbs form of the equilibrium states. 

The Jacobian matrix of the transformation relating the two parametrizations of T) is the entropic matrix 


d^S{R} 
dRa dRp 


dJ^ . n dHnZ{J} 
dRp' dJ°^dJP 


dRa 


(a,/3 ^ 0), 


(27) 


which is the (negative) matrix of second derivatives of S{R}. 

As a remark, we note that a metric ds^ can be defined in the full space of density operators D by ds^ = —d^S{D} 
[^ . The quantity —SRaS^^^SRp can indeed be interpreted as the square of the distance, within the trial subset of 
density operators, between the state V parametrized by {i?} and the state V + SV parametrized by {i? + 6R}. The 
matrix —S can thus be regarded as a metric tensor, and the relation 6J°‘ = —(a, /3 ^ 0) as a correspondence 
between covariant and contravariant coordinates. 


B. Symbols and images 

We shall have to deal with quantities of the form TrQI), where V will be some element of the Lie group and Q 
some operator in the full Hilbert space Jtf, not necessarily belonging to the Lie algebra. In order to take care of such 
operators it appears convenient to represent them by means of two useful tools. 

Let us first introduce the symbol q{R} oi Q, a scalar which depends both on the operator Q (in the Schrddinger 
picture) and on a normalized running element T) of the Lie group. This symbol is defined by 

q{R} = TrQV, (28) 

as function of the parameters Ra (with a ^ Q) that characterize V. If Q belongs to the Lie algebra, this function is 
linear since Ra = Tr Mq, V itself is the symbol of Mq,. Otherwise (/{i?} is non-linear. If V were a density operator, a 
symbol would be an expectation value but here T) is an arbitrary normalized element of the Lie group, not necessarily 
hermitian. The symbol q{R} of the operator Q may be viewed as a generalization, for any Lie group and for mixed 
states, of the expectation value of Q in a coherent state [30l| . this value being regarded as a function of the parameters 
that characterize this state. 

Let us then introduce a second object associated with an operator Q, its image Q{i?}, an element of the Lie algebra 
depending again both on Q and on the running operator R. The image Q{i?} of Q is constructed by requiring that 
both operators Q{R} and Q should be equivalent, in the sense that 

TvQ{R}V = TTQV = q{R}, (29) 

and that 

TtQ{R}SV = TtQSV (30) 

for any infinitesimal variation SR around R within the Lie group, in particular for SR oc MqIJ, or SR oc R Mq,, or 
SR oc dR/dRa- Let us show that these conditions are sufficient to determine uniquely the image Q{i?} associated with 
a given Q. Since it must belong to the Lie algebra, Q{i?} is parametrized by a set of coordinates Q°‘{R} according to 

Q{i?} = Q“{i?} Mq . (31) 

For a ^ 0, these coordinates Q“{i?} are determined by inserting (I5T]) into (1501) . by taking SR oc [dR/dRp)SRp and 
by using the definition (l28l) of the symbol of Q, which yields 

dq{R} 
dRa 




(a ^ 0). 


(32) 
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The coordinate Q°{i?} is obtained from (l29l) . Altogether, the image Q{R} of Q defined by their equivalence (I29t . (p0l) 
is related to the symbol q{R} of Q through 

Q{i?} = g{i?}Mo + (M„-i?„Mo)^|^ {a ^ 0). (33) 

C'./Xq' 

An operator Q belonging to the Lie algebra coincides with its image. Otherwise, its coordinates depend on the 
Ac’s, so that Q{R} is an effective operator simulating Q in the Lie algebra, for a state close to V. The dependence 
of on {R} arises from the occurrence of V in the equivalence relation (23) , (EOj) ■ 


C. The commutation matrix C and the entropic matrix S 

We shall have to handle products of two operators of the Lie algebra. When solving the dynamical equations we 
shall in particular encounter commutators [Mq,, M^g]. Their symbol defines the commutation matrix 

C<,p{R} = j^Tr[U^,Mp]V = Tl^R,, (34) 

expressed in terms of the structure constants T))^ (for 7 = 0 , we have i?o = !)• The matrix C will play a crucial role. 
Using lnT> = the vanishing of Tr[M^, In 2?] 2? = 0 is expressed by 

Rsr = 0 . (35) 

Taking the derivatives of this identity with respect to Ra and using the relation (1271) between the two parametrizations 
{J} and {27} of 27, one obtains for the product CS the relations 

= (a^O) (36) 

r = -Cp^n^^ Rs. 

These equations help to show that the automorphism of the Lie algebra engendered by the element 27^ of the group 
can be expressed in the two equivalent forms 

M„-i?„^27-^(M„-2?„)27^= (e*'^C®^)J(M^-i?^). (37) 

(This is proved by evaluating the derivative with respect to A of the left-hand-side, using (1361) then integrating from 
0 to A.) 

The property (1571) will be exploited later on ISec. lVEl and Appendix A). We use it here to find the expression of 
the symbol of the product M^, or equivalently of the correlation Tr — Ra Rp = Tr(Ma — Ra)Mpf> (with 

Ra = Tr Ma27). To this aim, we start from the expression (1771) of S“^, and evaluate explicitly therein the derivatives 
with respect to {J}: 

_d^\nZ{J} _ d Tre'^^'^^M;? 

)ap- dJaQjP “ Tre'^'"'^^ 

= Tr / d\V^-^{Ua-Ra)'D^Mp, (38) 

Jo 

where we made use of the first-order expansion in the shift {(5J} of 

exp{J'^M.y + 6rM..^)^V + 6J'^ [ dXV^-^M.yV^, (39) 

Jo 

with 27 = exp J''' M.^. We recognize in the r.h.s. of (l38l) the Kubo correlation of and in the normalized state 
27. By means of (I37p the integration over A can be performed explicitly in (1551) . which yields 

- (S“^)c /3 =Tr [ dX {M^-R^)Mpf> 

Jo 

/ zncs_ 7 

= A' 


( 40 ) 
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Hence the ordinary correlations between operators of the Lie algebra in any element T) of the Lie group are found to 
be given by 


Tr Ma MpT) — Ra Rp = I — 


gifiCS _I' 


5-1 


(41) 


D. The variational formalism in the restricted space 


With the above tools in hand, it is possible to implement the Lie-group form (I22p of the trial operators into the 
variational expression (IT^ of '^{A^'D} by taking as variables the parameters {L}, {^}, Z that characterize A and V 
at the times r or t. These operators A and V appear within T)} through products AV and DA. Such products 
belong to the Lie group and are characterized by the parameters R^^ = Tr Vh^AT)AT>R^'^ = Tr MaP^/Trl?^, 
= TtAV, which should be expressed in terms of our basic variables {L}, {i?}, Z by relying on the group 

properties. 

The initial state operator K = InZl, the Hamiltonian H and the observables Q® enter V} through traces 

of the form Ty QAT>/Tr AT> and Ti QV A/TyT) A, where Q stands for K, H or Q®. We are thus led to introduce, for 
any operator V parametrized by {i?}, the symbols 

k{R} = TYK'b, h{R) = TYH'b, = TrQ®^ (42) 

of K,H and Q®. These symbols occur within '^{A^'D} for values of {i?} equal to or {RP-^}. The variational 

expression '^{A,'D}, when specialized to a Lie group, then takes the form 


V} = In Tr A{t = M) V{t = /3) 
rP 

— dr 


(43) 


/ dlnZ-^^ dRa dlYiZ i r 


which should be regarded as a functional of the original trial parameters {L},{i?} taken at times r and t, and 
Z{t) = TyV{t). 

The stationarity eonditions (I15m8|) , obtained by functional derivation with respect to these parameters, now read 


dV{T) 

dr 

dAp) 

dr 

dAit) 


+ K{R^^)V{t) =0, 

- A{t) = 0 , 

1 


(0 < r < /3) 

(0 < T < /3) 


dt At)] 

+ tYl^:iAAt)Qj{R'^A = t^y it>ti) 

3 

\ AA H{R-^A - H{i?^-^}T'(t)l 
at in ■' 

-iJ2^jit)Qj{R^AAt) = o, it>ti). 


(44) 

(45) 

(46) 

(47) 


These equations involve the images K{i?}, H{i?} and Qj{7?} issued, according to the general relation (1551) . from the 
derivatives of the corresponding symbols k{R}, h{R} and qj{R} defined by (1^ . The stationarity conditions (l44l) - (l47l) 
should be solved with the boundary conditions ([T^ and m- When they are satisfied, is constant in r and t, 
and the optimal estimate for the generating functional ipi^} reduces to the first term InZ-^® of (|4^ . in agreement 
with (fTOl) . 

For most Lie groups, solving the coupled equations (|44II47I) is hindered by the need of expressing explicitly the 
parameters of AV and DA in terms of those {{L}, {i?} and Z) of A and V. However, we are 

interested in the first terms of the expansion of ^’{0 in powers of the sources ^j{t). Accordingly we shall only need 
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to express, as functions of the parameters {R} and Z of an arbitrary element T) of the Lie group, the following 
ingredients: (i) the symbols of the operators K, H and Q®, and (ii) the entropy function (1^ which allows us to 
relate the sets {i?} and {J}. This is feasible for many Lie groups. [In the example of the algebra of single-fermion 
operators this is achieved by Wick’s theorem.] Thus, explicit solutions of the equations of motion will be found 
at the first few orders in the sources. 


IV. ZEROTH AND FIRST ORDERS 
A. Thermodynamic quantities 

At zeroth-order in the sources {^}, the quantity of interest is the partition function Tre“^^, or the ’’generalized 
free energy” 

F=-/3-MnTre-'^^ = -/3-V{C = 0} (48) 

(/cb = 1,7’ = fi~^) which reduces to the standard free energy for K = H, or to the grand potential ior K = H — ^N. 
It is variationally approximated by — I3~^ InTr (t = (t = /3), where the upper index 

denotes the order in the sources {^}. 

For {^} = 0, the stationarity condition (l46|) yields = I ior t > ti, hence, from Ea. ifMl) . = j3) = I 

and F ~ — /3“^ lnTr2?(°^(T = /3). Thus = /3) = appears as an approximation, variationally suited to the 

evaluation of thermodynamic quantities, of the exact state D = e~ ^ ^. 

To obtain (r = /3) we have to solve the first two stationarity conditions (I44II45I) with (0) = (/3) = I. 

We make the Ansatz 

V{t)A{t) = , (49) 

where 2?^°) = (r = /3) is a constant operator, still to be determined and characterized by its parameters Ra^ = 
Tr for a ^ {) and = TrT>^^\ The image is then the constant operator so that we can 

solve Eas. dTil) and (1451) in the form A*'°^(r) = [T>(o)j(^-r)//3^ where is determined by the 

equation 

lnT>(°) =-/3K{i?(°^}. (50) 

The operator equation (1501) provides = /3) = More explicitly, in the basis {M} of the Lie algebra, the 

components a ^ 0 of (l50l) read = — (3IC°^{R^^'>} in terms of the coordinates of ln2?(°) = and the 

coordinates of [defined as in (EIllSSl)]. This yields the self-consistent equations 

dS{R^^^ _ ^dk{R(°'>} 

which determine the parameters {2?^°^} of 'D^^\ The component a = 0 of (ISOl) yields — 

Ra^ where Eg. (1001) has been used for K. Together with (1001) and (15T1) . this equation provides for the 

sought free energy F the alternative expressions 

F ~-/3-MnTre-^'^{^'°’> =-/3“^ , (52) 

where we defined the free-energy function f{R\ through 

f{R} = k{R}-TS{R}. (53) 

The above self-consistent equations determine a hermitian operator 2?^°^ which can be interpreted as an approx¬ 
imation for the exact density operator D (whereas the trial operator 2? occurring in the presence of sources is not 
hermitian). The general relation (1051) entails 


=0, (54) 

an equation which is equivalent to Tr M.y[ln2?(°\ = 0. [In the Hartree-Fock example, ((50)) expresses the com¬ 

mutation of the single-particle density matrix with the effective hamiltonian matrix.] 
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Our variational principle provides solutions for il}{^ = = — P F that are not maxima but only stationary values of 

2?}; it relies on the Bloch equation rather than on the maximization of the entropy under constraints. However, 
it turns out that the above result (EH coincides with the outcome of the standard maximum entropy procedure. 
Indeed the latter amounts to minimizing the left-hand side of the Bogoliubov inequality 

TtKV + > -/3-MnTre"'^^ = F (55) 


with respect to the normalized trial density operator T>. When T> in (l55]l is restricted to the Lie group, the left-hand 
side reduces to f{R}, where {i?} is the set characterizing V. The best estimate for F is thereby the minimum of 
the free-energy function f{R}, which requires that the equations (ISTI) are satisfied. The equivalence between the 
two variational approaches provides a criterion for selecting the best solution of the self-consistent stationarity 

conditions (1511) when they have several solutions, namely, the one for which f{R^^^} is the absolute minimum of f{R}- 
The standard relations S = —dF/dT ~ — Tr(D(°) ln(D(°) = i^) — ^ + TS ~ are satisfied 

as usual, so that the approximation is thermodynamically consistent. Thermodynamic coefficients are obtained by 
derivation, which introduces the matrix (a, /I ^ 0) 


dR^P dR^P dRf ' 

In particular, the heat capacity is found for K = FI as 

C 7 ~ k:<^{r^°'>}. 


(56) 


(57) 


The positivity of the stability matrixW at the minimum of f{R}, entailed by the inequality (I55L will play an important 
role below. [For the fermionic single-particle Lie algebra we recover the thermal HF approximation, either under 
the form (1501) or through the minimization of f{R}.] 


B. Expectation values 


Expectation values are obtained by expanding the generating functional i/'IC} fi^'st order in the sources ij(t) 
[Eq.([5])]. These sources occur both directly in '^{A^V}, as exhibited by the last term of (|i5)) . and indirectly through 
the values of A and V at the stationarity point where V’lC} = R}- At this point, however, the partial derivatives 

of 4'{^,27} with respect to A and 27 vanish; one is left with the explicit derivative d'^{A^^\'D^^^}/d^j{t) taken at 
the zeroth-order point {^(°(,27(°)}, so that 


{Qj)t 


IdiPiO TrQf27(0)(t)yl(0)(t) 

“ Tr^(o)(t)27(o)(t) 


(58) 


involves only the zeroth-order approximations ^(°-*(<) and 27(°)(t). 

We have seen that A^^^t) = I ioi t > fi, and that 27(°)(ti), equal to 27(°)(r = P) = is given self-consistently 
by (|5ni) . Hence, the static expectation value {Qj)^. in the state D is variationally expressed by 


Tr Qf D = - Tr Qf , (59) 

that is, by the symbol (l28l) of the observable Q® of interest. 

For dynamical problems, the expectation value (I58|) involves 27(°)(t), provided by the last stationarity condition 
EH taken for ^j{t) = 0 and for ^(°(’(t) = /. Thus, 27(°)(t) is determined by the zeroth order self-consistent equation 


dV^°\t) 

dt 


H{R^°\t)}, V^°\t) 


(60) 


with the initial condition 'b^^\t{) = [The norm Z^'^'>{t) is constant and equal to Z^^\] In (|60)) the operator 

H{i?(°(*(t)} is the image taken for Ra\t) = Tr Ma27(°)(t) of the Hamiltonian 22. The time-dependent expectation 
values (155)) are therefore found as 


{Q,)t = TvQfiptOD ^ TrQSp(o)(i) = q^{R(o){t)} . 


( 61 ) 
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As noted at the end of Sec. uni the variational equation (l60l) for the Lagrange multiplier comes out as an 

approximation for the Liouville-von Neumann equation (1201) . 

In coordinate form, the equations of motion for the variables Ra \t) = (Mq,)* that parametrize are found 

from (l60ll as 


dR^a\t) 

dt 




dh{R(°\t)} 

dR^it) 


(62) 


with the initial conditions {R^^\ti)} = [For the single-fermion Lie algebra, we recover the time-dependent 

Hartree-Fock (TDHF) approximation for the single-particle density matrix R!^^{t), with the static HF solution as 
initial condition. Indeed, the multiplication in (1621) by ^lap{R} produces for this algebra a commutation with R^il{t). 
The usual single-particle effective HF Hamiltonian comes out as the image The current use of TDHF to 

evaluate expectation values is thus given a variational status.] 

Equations (1621) can be rewritten in the alternative form 


d^ 

dt 




(63) 


where the matrix 


LJ{i?(o)(t)}=rJ^7^^{i?('’)(t)} (64) 

plays the role of an effective Liouvillian acting in the Lie algebra; the quantities are the coordinates = 

dh{R}/dRi 3 [Eg. (15^ ] of the image H{i?} of H. 

The dynamics of the density operator {t) takes therefore the classical form (l63ll in terms of the scalar variables 
Ra^ parametrizing This classical structure will be analyzed in Sec. B Non-linearity occurs through the 

restriction of the trial space: the image H{i?} depends on the variables {R^^^t)} when H does not belong to the Lie 
algebra. 


V. CORRELATION FUNCTIONS 

Variational approximations for the two-time correlation functions Cjk(t',t") are provided by the second-order term 
of the expansion ([8]) of the generating functional in powers of the sources However, we have seen that, 

owing to the stationarity of 4/, zeroth order was sufficient to determine the variational approximation for expectation 
values. Likewise, it is sufficient here to expand up to first order the first derivative of the generating functional 'i/'ICI 
according to 


15 ^ 

i 


{Qk)t" “f ^ 



3 


\d^{A,V} TTQlV{t")A{t") 
~i d^k{t") ~ TrV{t")A{t'') 


(65) 


where the stationarity of '^{A^'D} with respect of A and V has been used. We can thus obtain Cjk{t',t") from the 
first-order contribution to the r.h.s. of (p5t . and only A^^'^ and have to be determined from the 

variational equations (I44M7I) . 


A. The approximate backward Heisenberg equation 

One building block for correlation functions will be the quantity Q^(t",t) defined for t” > f by expanding the trial 
operator A{t), which enters (l65l) . up to first order as 

poo 

Ait)^A^°\t)+A^^Ht) = I + ^ / dt"Y,UnQkit''A)- 

Jt 


(66) 









14 


Comparison with the expansion of the exact generating operator A{t) defined by (|4]) shows that simulates 

the Heisenberg observable Since A{t) belongs to the Lie group, belongs to the Lie algebra and 

can be expressed as 

The stationarity condition (H51) with respect to A{t), expanded up to first order, determines the coordinates 
For a ^ 0 these obey for t < t" the equations 

= (a,/3^0), (67) 

which appear as the reduction in the Lie algebra of the backward Heisenberg equation (11011 for Q^it",t)- The matrix 
L is the effective Liouvillian (IMl) , C is the commutation matrix defined by (IM)) and H is the second-derivative matrix 
of the symbol h{R} of the Hamiltonian iL, 




d^h{R} 

dRa dRp 


( 68 ) 


These three matrices are functions of {i?}; in (1571) they are taken at the point Ra = Ra\t) determined by the 
zeroth-order equations (l63l) . Equations (l67l) should be solved backward in time from the final boundary condition 


QT{t"X) = Qt{R^^\t")} 


dqk{R^^\t")} 

dR^°\t") 


(69) 


where Qk{R}^a = Qk{R} is the image of the observable Qk while qk{R} is its symbol. 

The exact Ehrenfest equation for Tr MaD{t) [issued from the Liouville-von Neumann Ea. (l20l) ] and the exact Heisen¬ 
berg equation for Q^{t",t") involve the same Liouvillian. However, in the variational treatment, the corresponding 
approximate Ehrenfest equation (I63|) and the approximate Heisenberg equation (16711 differ; the latter contains the 
corrective term CH in addition to the effective Liouvillian L. This difference arises because the variational equations 
for the set Ra\t) occur at zeroth order in the sources whereas those for the set Qf“(t",t) occur at hrst order. 

From the stationarity condition (l46ll for A{t) one also Hnds the component Q™{t”,t) of . This gives an 
alternative expression for the time-dependent expectation value of Qk at a time t, 

{Qk)t qk{R^°\t)} = = TrQ^{tA")V^°\t") (70) 


(including a = 0 with = 1), which holds for any intermediate time t", and thus interpolates the Schrddinger 

picture for t" = t [Ea. dbTll ] and the Heisenberg picture for t" = ti, as does the exact expression 

{Qk)t = TrQf{tA")[U{t’'Ai)DU\t",U)]. (71) 


B. Bypassing T>^^l{t) for real times t 

The correlation functions that we want to determine through (l65)l depend on the combination 'D^^\t") + 
A^^\t"), since A^^'^ = X. The time-dependence of A^^\t") has been expressed through Eos. (1551) and (1^ . 
We still need the other ingredient, the first-order operator T>^^)(t") whose evolution is not simple. However, using the 
coupled equations of motion for V{t) and A{t), we will show below that the time t” occurring in (1551) can be shifted 
down to ti. 

From Eas. (l451) and (l47l) . one can derive for the products 'D{t)A{t) and A{t)'D{t) the uncoupled equations 

^ RA] + [QAR'^^}, RA] , (72) 

3 

^ AV] . (73) 


These equations cannot be fully solved because the boundary conditions on A and V occur at different times [Eg. dTSll ]. 
Nevertheless, Ea. (l72l) will happen to be sufficient for our purpose. 
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In parallel with Ea. (l66l) . we parametrize at first order the product V{t)A{t) by according to 

poo 

= R^^\t)+z dt'J2^,{t')R%\t,t'). 


( 74 ) 


We have in (1741) denoted as SV^t) the first-order variation of 27^/TrX>^ around The quantity Cjk{t',t") that 

we wish to evaluate now satisfies from (1551) and (174)) the relation 


poo _ 

+ t dt'Y] Cjk{t', t") = Tr QS sv{t") . 

Jti 


(75) 


Since S'D{t") is a variation around within the Lie group, the image property (l30l) allows us to replace the 

operator Q® by its image Qfc{i?(°i(t")} = {i?(°i(t")}M^. Comparing Eos. (1741) and (1751) . we can get rid of the 

sources so that 




(76) 


Remember that the coordinates of the image of are related to the symbol qk{R^^\t'')} = 

TrQSp(o)(i//) by 


dR^it") 


(77) 


The time-dependence of R^j^J (t, t') is found by expanding Ea. dT^ at first order. Noting that TiVA does not depend 
on time, one obtains 


dR^^Ht t') 

’ = (L + CH)^^ (t, t') + thC^^ Qj{-R^°^(t)} Ht - t'), 


(78) 


where again L, C and MI depend on time through {i?i°i(f)}. The same kernel L-bCH as in the approximate backward 
Heisenberg equation (1671) is recovered; it will also be encountered in the context of small deviations ISec. lIX Bll . [For 
fermionic systems, L -|- CMI is the time-dependent RPA matrix issued from the kernel L of the TDHF equation (1541) .] 
As a consequence of the duality between the kernels of Eas. dHTl) and (1781) . we obtain the identity 


t)] = -zh Q]{R^°Ht)} c,p 6{t - t'). 


(79) 


Since Q^^(t",t) vanishes for t” < t, the r.h.s. of (174)) disappears if t" < t'. One can therefore evaluate Cjk{t',t”) for 
t' > t" from ()751) by noting that the product Rf'Jit, t') Q^^(f", t) does not depend on t for t' > t” > t > U. Using the 
boundary condition Q”'^(f",t") = Qi{R^°\t'')}, we then shift t" down to t; in (174)) . which yields 

The continuity condition T>{ti)A{t{) = 'D{I3)A{P) entails that this correlation function (1541) is equal to 

= {t'>n ( 81 ) 

in terms of the boundary value at r = /3 of {i?j^^(T, t')}. One needs therefore to determine, only in the interval 
0 < T < /3, the set {i?j-^^(T, t')} defined through the first-order parametrization of V{t)A{t) by 


Tr T’('r) A{t) 
TvV{t) A{', 




R^^ir) 


(82) 
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where i?® = Tr . 

The evaluation of at times t > ti has been bypassed, and it remains to solve the coupled equations (l4^ and 
(HSl) for and in the range 0 < t < ^, with the boundary conditions = 0) = 0 and 


= (3) 


= i 


/ CO 


(83) 


issued from m and dMl), respectively. For t” > t', the symmetry Cjk{t'A") = Ckj{t”,t') can be checked by means 
of the r.h.s. of Eg. (1751) . 


C. Two-time correlation functions 

The linear structure of the boundary condition (I83p for together with the boundary condition = 

0 ) = 0 for entail the occurrence of an overall factor in A^^\t) and hence in R‘d^{T,t'), 

and finally in the correlation function Cjk{t', t"). We had already acknowledged the explicit dependence of Cjk{t', t") 
on Therefore, for any trial Lie group, the two-time correlation functions have the general form (a,/? ^ 0) 

Beside the approximate Heisenberg observables given by the dynamical equations (1671) and the boundary condition 
dMl), a matrix appears, which depends only on the zeroth and first order contributions to T’(/3) and A{f3). 

Noticeably, the expression (1841) has the same factorized structure as the exact formula dZ]). In spite of the coupling 
between the stationarity conditions in the sectors t > U and 0 < r < /3, Eg. (15^ [as Eq. dZl] displays separately two 
types of ingredients: The result B of the optimization of the initial state is disentangled from that of the dynamics, 
embedded in the approximate Lie-algebra Heisenberg observables Q^{t'Ai) and Q^{t" Ai)- 


D. The correlation matrix B 


The last task consists in determining explicitly the correlation matrix Bq,^. We have noted above that, through 
the boundary condition (1551) . the factor Qfj°‘{t'Ai) occurs in i?j^^(r,t'). By introducing new (a-indexed) operators 
A^a\T) and T>a^(r), again determined by the coupled equations (l44l) and (1^ but with the boundary conditions 
Aa\l3) — and 'Da\l3) = 0, one can explicitly factorize R^-^p{TA') as 

i?;^^(r,0 = Q?“(i',ii)i?S(r). (85) 

The quantities R^^^ (r) defined by (|57|) and (1551) will then be given by 

i?«(r) = Tr[2?W(r)^W(r) +pW(r)^(o)(T)](M^ - , (86) 

and the correlation matrix Bq,^ will be found as Bq,^ = = /3), that is. 




ivpi.”(/i)(Mg - R'; 


(0) 


(87) 


It remains to determine the coupled operators Va^ (r) and Aa^ (r). The solution is explicitly worked out in Appendix 
A.l. It had already been given in special cases, for fermionic systems extended BCS theory [13, (jA quantum field 

M- 


The result thus found for the correlation matrix B is 

ihC¥ 


F“ 


1 — exp(—i CF) 
iinCFcoth(ii;i/3CF) 


(88) 


¥-^ + -ihC. 






17 


It involves the commutation matrix C defined by (p^ and the positive matrix F = K — T § of second derivatives (l56l) 
of the trial free energy. Both are taken for the parameters that were determined at zeroth order. For vanishing 

eigenvalues of iCF, the coefficient of F“^ is meant as [For the single-particle fermionic Lie group, the matrix 

iCF is the static RPA kernel.] 

By letting t' — 0 — t" = ti and Qf = Maj in (l84)) . one identifies as the variational approximation for 

the correlations, in the initial state D, of the operators Mq, that span the Lie algebra: 

TrM„M,3Z)-TrM„Z)TrM;3L>=iB„;3, {a, ^ ^ 0). (89) 

The variational expressions ([84]), (I88|) [together with Eas. (l67)) ] for the correlation functions issued from the use of 
a Lie-group trial space are the most important outcomes of the present variational approach. We comment below the 
features of this expression, and work out its consequences in the forthcoming sections. 


E. Status of the result for static correlations; Kubo correlations 


The optimization of thermodynamic quantities and expectation values (^Sec. lIVI) has resulted in a mean-field type 
of approximation, with the mere replacement of the exact state D by the zeroth-order contribution 2?^°^ to the trial 
object V, as in Tr MaD ~ Tr = Ra^. However, the optimized approximation (I88L(I5^ that we found for the 

correlations Map of Mq, and does not follow from such a simple replacement. Evaluated in the state by the 
formula m, such correlations, instead of (1551) . would be given by 


TTMaMpV^°^ - Rf'> 


Vexp[ifiCS] - a 


(90) 


where S and C are evaluated from (1571) and (IMl) for {i?} = Contrary to B, the naive expression o is not 

variationally optimized. While the first term of the variational expression of B provides a contribution equal 
to ([501) . its second term, arising from introduces a correction which substitutes S — /3IK = —/3F to § . The 

matrix K takes into account effects coming from the part of the operator K that lies outside the Lie algebra. [When 
{M} is the fermionic single-particle algebra, the left-hand-side of (TOl) is the Fock term from Tr (aj)aj,)(aJ.aT-)27^°^ 
since R^^ is the Hartree term. The full matrix B involves an RPA kernel, in which the matrix K is the effective 
two-body interaction.] 

Although the approximations for expectation values and correlations functions stem from the same variational 
principle, they appear intrinsically different. Not only B cannot be expressed from 'D^^\ but moreover there is no 
density operator approximating D in the original Hilbert space that would produce the optimized correlations B 
in the same form as (l89l) . While can be interpreted as a state, the trial operator V has no perturbative status 
and is just a calculational tool involving the sources. In the expression (1571) of B, the operator 'Da \ which depends 
on Mq, is not a correction to 2?^°). 

Nevertheless, we will show in Sec. lVIll that, through a mapping of the original Hilbert space into a new space 
and of the Lie algebra {M} into a reduced algebra {M}j the matrix elements of B can be interpreted as exact 
correlations between the operators {M} in an effective state The quantities will also appear as exact 

expectation values of {M} in 22- The variational results, for both expectation values and correlations, will thus be 
unified through a modification of the Lie algebra. 

One may also wonder about the origin of the Bose-like factor exhibited, for arbitrary Lie groups, by the expression 
(155)) of B. A clue will be given in Secs. lVIll and IVlID where an algebra of Bose operators arises from the mapped Lie 
algebra {M.}- 

An alternative understanding of the structure of the matrix B can be reached by relating it to the matrix of Kubo 
correlations, even though the latter have less direct physical relevance than ordinary correlations. Kubo correlations 
between the operators Mq and M^ are defined for the exact state D by 


1 /■p _ _ _ 

Tr- / dre^^ MaC-'^^ MpD-TTMaDTrMpD. (91) 

P 7o 

We have seen in Sec. lHI Cl that, if the state D were bluntly replaced in (1971) by the element 2?^°) of the Lie group, 
Kubo correlations would be directly related by (1551) to the matrix §. We show in Appendix A.2 that a variational 
approximation B^ for the Kubo correlations (1911) reads 

B^ = ^(IF"^) = (/3K-S)-^ 


( 92 ) 
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The naive approximation (l38)) . namely ~ obtained by replacing D by in (EU, is thus variationally 

corrected by inclusion in (IMl) of the term (3 K. Likewise, the naive approximation (TOl) for the ordinary correlations 
in the state D results from the variational expression (|Mll for B by omission of K within F = K — T S. 

Once B^ has been obtained in the form (IM)) . it is possible to recover from it the expression (1551) for B. In the 
special case K = 0, for which /IF = —this was achieved in Sec. lIII Cl We saw there that the factor — 

entering the Kubo correlation (lOTl) is expressed as (M^ — i?^) through the automorphism 

dSZl) of the Lie algebra. Integration of between 0 and /3 yields 


,ihcs 


-I. 


ihCS 


in agreement with (I4II1 . In the general case K ^ 0, it is shown likewise in Appendix B that 


(93) 


B(B^)”^ 


ihpC¥ 

1 — exp(—i h/3C¥) ' 


( 94 ) 


This ratio stems therefore from a property inherent to the Lie group underlying our variational approach, namely the 
exponentional form of the automorphism (1371) . The Bose-like factor that enters B arises from this property, and from 
the simplicity of B^. [The quasi-boson structure of the static thermal RPA for fermions [la appears as a special case, 
see Sec. lVIIll ] 

Thermodynamic quantities and expectation values were determined in Sec. lIVI they depend on the mean-field images 
of K and H within the Lie algebra, namely, which self-consistently (Sec. lIV A)) determines and H which 

governs ISec. lIVBl) the time-dependence of More elaborate ingredients are required for the evaluation of 

static and dynamic correlation functions: The matrix IK enters B through F while the matrix HI (together with the 
mean-field Liouvillian L) governs the time-dependence of the approximate Heisenberg observables 

In the rest of this article we will analyse the properties of the approximate correlation functions [expressed by 
Eas. (l84|) . (|88ll and (|67|) ] and of the dynamics [expressed by Eqs.(|62|)]. 


VI. PROPERTIES OF APPROXIMATE CORRELATION FUNCTIONS AND FLUCTUATIONS 

In this section we review some consequences of the variational expressions found above for the correlation functions, 
encompassing special cases and conservation laws. 


A. Time-dependent correlation functions in an equilibrium initial state 

We first consider the special case of an equilibrium initial state D = exp(—/lA) for which the operator K is equal 
to the Hamiltonian H, plus possibly some constants of motion such as the particle number for a grand-canonical 
equilibrium. One can then generate simply the dynamics by using K instead of H in the backward Heisenberg 
equation (nnD. 

The exact expectation values {Qj)t ^^en do not evolve. Their approximation (ISTl) depends on time through {i?(°)(<)} 
governed by (|63ll . In the effective Liouvillian L defined by (l64)l . the coordinates of the image of the Hamiltonian 
are replaced by those of the image K of K, that is, = dk{R}/dRj. Using the self-consistent equilibrium 

condition (1501) for and the identity (|501) for the product C§, we find 

L^“{i?(o)} = = -Tr^“ = . (95) 

The identity (150)) then entails, as expected, that Ra\t) = Ra \ and hence that the approximation {Qj)^ does not 
depend on time. 

For two observables Qj and Qk that do not commute with AT, the exact two-time correlation function Cjk(t',t'') 
defined by o depends only on the time difference t' — t". This property is not obvious for the approximation (I84L 
which involves two factors depending separately on t' and t". To elucidate this point, let us solve for HI = IK the 
approximate backward Heisenberg equation (1571) for the image Qf{t',t) = Qf{t' — t) with the boundary condition 
Q]^(0) = QjfR^^Y- This equation involves the kernel L -|- CH which, according to (I95L takes the simple form 
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C(]K — TS) = CF. Hence, the Heisenberg equation (l67l) simplifies into ,t)/dt = — (CF)^, where C 

and F are evaluated for it is readily solved as 






(with a,/3 ^ 0) in terms of the boundary condition . 

In order to evaluate the correlation function (IM)) . use is also made of 


Qr(i",ti) = 






and of the relation 

F“1 g-®'C(t"-ti) _ g-CF(t"-ti) ]p-l ^ 

The explicit dependence on time of the correlation functions is then given for t' > t" by 




,CF(t'-t") 


ihC¥ 


1 — exp(—z ?i/3CF) 


F- 




(96) 


(97) 


(98) 


(99) 


The identity of H and K has led to the occurrence of the same matrix C F in the dynamical equation (|%1) and in 
the matrix B [Eg. dSSl) ] that accounts for the correlations in the initial state. As a consequence, only t' — t" appears 
in HMD, as it should. This property would not have been satisfied if the correlation matrix B were naively evaluated, 
as in Eq. o, by replacing in (IMl) the exact state D by 

For j = fc, Cjj{t' — t") provides the variational approximation for the autocorrelation function of the observable Qj 
in the equilibrium state D oc exp(—/SAT). 


B. Other special cases 

1. Commutators and linear responses 

The antisymmetric part of B, namely ih'C/2, is simple and depends only on the zeroth order in the sources, not on 
F. As a consequence the approximation Bq,^ — B^jq, for the expectation value Tr [Mq,, M/gjU = of the 

commutator [Mq,, M/j] is obtained as ihCa/ 3 {R^°^, a property in agreement with the expectation value R^'^ of 
found at first order in the sources. 

Linear responses are expectation values of commutators, and therefore involve only this antisymmetric part of 
Cjk{t',t"). Alternatively, they can be evaluated directly from the variational expression 'L{A,27} by including a 
time-dependent perturbation in the Hamiltonian. Then, the response of Qj to a perturbation Qk appears as an 
expectation value, and is hence directly obtained at first order in the sources, without the occurrence of K which 
enters the symmetric part of Cjk{t',t”). Both approaches yield 

xjkit\t")= ii/zh)0{t' - 

^ e{t' -1") Qf“(t',t)CQ4i?(°)(i)} , (100) 

where t is an arbitrary time in the interval U <t < t" and 9 the usual step function. In particular, by letting t = ti in 
(I100|) . the responses are variationally expressed in the Heisenberg picture in terms of the matrix and of the 

Heisenberg observables Q^{t',t{) and Q})(t",ti) given by the approximate backward Heisenberg equations (1571) . 

For H = K, D !x exp{—l3K), the response (llOOl) in an equilibrium state depends only on the time difference and 
takes the form 


Xjk{t',t")= {l/ih)9{t' -nTiD 




Qj ^ 


— i H {t' —t" 


Qk 


~ e{t' - t") d ■ 


0.0 


( 101 ) 
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2. Static correlations; classical limit 


Static correlations between observables Qj and Qk in the state D oc e 
t' — 0 = t" = ti in (15^ , which yields 


are variationally obtained by letting 


Cjkiti + 0,ti) = ■ (102) 

The matrix B of correlations between the operators {M} is here saturated by the coordinates = 

dqj{R^°^}/dRa^ of the images in the Lie algebra of the considered observables. 

For an Abelian algebra, the matrix C vanishes and the ratio (iMl) reduces to unity. The ordinary and Kubo 
correlations are identical, and the matrix B simplifies into = (/3F)“^. In the high-temperature limit /3 —>■ 0 and in 
the classical limit fi —>■ 0, the ratio (|M)l also tends to I, so that 

B^ (/3F)-i. 

If Qj and Qk are commuting conserved observables, the occurrence of the commutation matrix C in this ratio also 
reduces it to I. 

For a Curie-Weiss model of interacting spins aj = ±1 at equilibrium, the Weiss mean-field expressions for the 
thermodynamic properties and the expectation values (cr_,) are recovered, while the Ornstein-Zernike approximation 
M for correlations is recovered from B = (/3 F) 

For the fermionic single-particle Lie algebra, the matrix B is the variational approximation for the correlations of 
the operators and 0 ^. 0 ,-: 


Tt {al^a„){alar)D - Tr aj^a^DTr a^arD ~ ■ (103) 

If D were replaced by an independent-particle state V, the expectation value Tr (a^ ai/)(al ar)2? would be given by 
Wick’s theorem and Tra^ 0,^22 Trala^^ would be the Hartree term so that the left-hand side of (11031) would reduce 
to the Fock term Tr(al 0x2?) Tr(a,yal2?). For a more general state 22, the expression (IMll of B involves the static RPA 
kernel 2CF, which takes into account not only the Fock term but also, through K, effects of the interactions present 
in K. 


3. Fluctuations 

The static fluctuation AQj of the observable Qj in the state D is variationally given by (11021) where k = j, that is, 

Aq 2 = Qf{2?(°)} . (104) 

This fluctuation AQj{t) evolves in time according to (1841) with k = j, t' = t" = t, that is, 

AQ2(i) , (105) 

which involves only the symmetric part of (1551) . For H = K, the fluctuation is time-independent. For H K, 
examples of time-dependences that are not properly accounted for by non-variational mean-field approximations are 
given at the end of Sec. IVI Al and in Sec. IVI Cl 


4- Initial state in the Lie group 

In case the exact density operator D belongs to the trial Lie group, the operator K belongs to the Lie algebra and 
coincides with its image, equals 22; the matrix K vanishes, F reduces to —TS and B to the trivial form (1901) . [For 
fermions, B reduces to the Fock terms.] The quantities ti) and Q^^(t", ti) are in this case the only ingredients, 

apart from 'D^^\ that enter Cjk{t',t"). 
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5. Zero-temperature limit 

The present formalism encompasses ground-state properties, found by letting /3 —>■ oo. This limit entails simplifica¬ 
tions. [For instance, in the fermionic case, the parameters of constitute a matrix satisfying [1?^°^]^ = 

The number of vanishing eigenvalues of the commutation matrix C increases. While S{R^^^} tends to 0 as T —?> 0, 
the quantities /3~^dS/dRa \ /3“^S and CF remain finite, due to the singularity of the von Neumann entropy (1^51) for 
vanishing eigenvalues of V. The resulting simplifications of the correlation matrix B will be exhibited below in the 
diagonalized form (jl35|l of B. 

A further simplification occurs if the initial state lies in the Lie group. For fermions possibly with pairing) the 
initial state is then a Slater determinant (or a BCS state). In this case, it has been shown (6| that one can by-pass the 
solution of the equations (l67)) for Qf°‘{t',ti), of the RPA type. The proof relies on the fact that these equations (l67l) 

for Q^°‘{t',ti) involve the same kernel as the dynamical equations (11651) for small deviations SRa\t), and that the 
latter equations can in practice be worked out by expansion of the simpler time-dependent mean-field equations (16211 
for Ra\t). Two-time correlation functions Cjk(t',t") and time-dependent fluctuations AQj{t) can thus be evaluated 
by running the existing TDHF (or TDHFB) codes alternatively forward and backward, with appropriate shifts in 
the boundary conditions. (For another derivation, see @,[110 .) This technique, variationally consistent, has been 
successfully applied [10 to describe, in nuclear systems, the fluctuations of single-particle observables which were 
severely underestimated by the conventional use of TDHF; for a review, see [HI- 


C. Images of Heisenberg operators; conservation laws 

It has already been noted that the equations of motion (1571) for the coordinates Qf°‘{t',t) of appear as 

a variational counterpart of the backward Heisenberg equation m- To be more precise, let us write the time 
dependence of the Lie algebra operator Q^(f', t) = Qj^“(t',t) M^. Using the relation (L-kCIH)^ = T^^'H'*{R^°^t)}-\- 

C,g.y{i?(*’^(t)}]HI'''“{i?(°)(t)} = d{<Cp.y'}V)!dRa \t) and the definition (155)) of images, one can rewrite the equations of 
motion for Q^“(t',t) (including a = 0) as 

dQ^it^ t') r 1 I 1 

— ^’ = Image of — [Q“(t',t), 7J] j = - — [QH(t',t), H] , (106) 

the image H of the Hamiltonian H being evaluated with respect to {R^^\t)}. While the equations (1571) were written 
in terms of the coordinates of QP, the introduction of images gives these equations a simple operator form: The 
right-hand side is simply the image of the r.h.s. of the exact backward Heisenberg equation m for time-dependent 
observables. The boundary condition Q^{t',t') = Q® is also the image of the boundary condition for the Heisenberg 
observable Qj{t',t)- 

Up to now, no specific assumptions have been made about the data. Let us now consider an observable Q® that 
belongs to the Lie algebra and commutes with H. This conservation law, together with the equation of motion ()106p . 
shows that the operator Q^{t',t) is constant and equal to Q®. Hence, the approximate expectation value {Qj)t, 
evaluated through (ITOl) for t' = t;, is constant. Less trivially, the fluctuation AQj{t), evaluated through (I105p . is also 
constant as it should. This property, which arises naturally through the use of the approximate backward Heisenberg 
equation, was not granted. [For instance, in the time-dependent Hartree-Fock approximation for fermions, fluctuations 
evaluated from through Wick’s theorem are not constant for a conserved single-particle observable.] 

Another conservation property holds for two observables Q® and of the Lie algebra such that Q® is conserved, 
[Q®, H] = 0, and that [Q®, H] = ihQj. This occurs for instance if Q® = X is a coordinate of the center of mass 
of the system, and Q® = P/m the associated velocity operator (m is the total mass). Then, Ea. ()106l) implies that, 
as the exact Heisenberg operators, the approximate ones satisfy P^(t',t) = P® and X^(t',t) = X® -f (t' — t)P®/TO 
(see Sec. 5.5 of i)- Hence, not only (X)^ and (P)^, but also the fluctuations and correlations of X and P produced 
by the general formula (l84l) have the proper time dependence. In particular, AP is constant and AX^(t) satisfies 
dAX'^(t)/dt = AP^/to^; we thus acknowledge that the present approximation for the fluctuations accounts for the 
exact spreading of the wave packet [whereas the TDHF approximation produces a constant width]. 
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VII. MAPPED LIE ALGEBRA AND MAPPED HILBERT SPACE 

The purpose of this section is to rewrite in a unified form the variational approximations (I59|) . (I6H) for expectation 
values and (l5^ . (|88p . (IMl) for correlation functions . To this aim, we will rely on a correspondence that associates with 
the Lie algebra {M} a simpler Lie algebra {M.}- It will turn out that the various expressions found above can be 
re-expressed as traces over a single effective density operator D, acting in a new mapped space rather than in the 
original Hilbert space 


A. Unifying the approximate expectation values and correlations in the initial state 


The optimization of the expectation value (M„) = TrllMa has yielded the approximation (MQ)g^pp = = 

TrD(o)M a', in contrast the optimization of (M„M^) = Tr_DMaM^ has yielded (MqM/3)^pp = + Ra^ that 

cannot be expressed as a trace over a density operator in the Hilbert space We wish to map the set {M} acting 
in Jif onto a new set {M} acting in a new space and to introduce in an effective density operator D so as to 
re-express our approximations in terms of D_. Namely, we wish the exact expectation values over D, in the mapped 
space denoted as (Mct)map and (MaM^)map, to coincide with the corresponding variational approximations in the 
original space denoted as (Ma),jpp and (MQM/3)app, according to 

(Mc)map = 2lMo D = (M„)^pp = , (107) 

(McM/3)map = TLhLMpD = (M„ M;3 )^pp = . (108) 

Going from the space ^ to will be a price to pay for expressing both expectation values and correlations in terms 
of a unique effective state 

The first step consists in replacing, in the original Lie structure [M^, \\Ap\= i M.y, the operator M.y on the right 

side by the c-number r!^'^ = Tr the expectation value of in the state of M’. This procedure associates 

with the original Lie algebra {M} in a reduced Lie algebra {M} characterized by the simpler commutation relations 

[Mc,Mp]=ifir20R^°^M.o = ifi‘Co.p. (109) 

(From now on we shall most often drop, as in the end of (jl09|) . the unit operator Mo-) 

Multiplication of any number of operators M ^. generates an enveloping algebra, the space of representation of which 
defines The structure of this space will be cleared up in Sec. lVHlI bv setting Cap into a canonical form. 

In order to satisfy the conditions (I107p on expectation values (Mc,)„^pp the sought effective density operator D 

should depend only on the differences Mq ~ Ra^■ As regards the conditions (I108|) . we remember that correlations 
are often generated in statistical mechanics by a probability distribution having the form of an exponential of the 
free energy regarded as a function of the running variables. For instance an energy distribution is the product 
e~^F{E) _ ^-p E+S(E) qJ |.]^g Boltzmann-Gibbs exponential e~^F by the level density, the exponential e^^F) Qf ^^.e 
entropy, which accounts for the discarded variables. Another example was provided in our approximation by the 
thermodynamic coefficients given by the second derivative (1551) of the free energy f{R} at {i?} = {R^^^}. Here, 
likewise, the reduction from {M} to {M} suggests to rely on an effective free-energy operator rather than on an 
effective Hamiltonian. We therefore replace, in the free-energy function f{R} = k{R} — TS{R} defined by (l53l) . the 
variables {i?} by the corresponding new operators {M}- The operators Mq fluctuate around Ra \ and this leads us 
to expand the operator /{M} as 

/{M} ~ f{R^°'>} + ^(Mq - A®) F“^ (M^ - <)) +..., (110) 

where the first-order term is absent thanks to the stationarity of f{R} at 

We thus guess that the distribution D governing the operators {M}; which is expected to yield the identities 
(I107l) - (|108l) . should have in the mapped space the exponential form 

D= Z=i(MQ-i?®)F“/5(M^-Af), (111) 

where £ behaves as a kind of free-energy operator. The inclusion of an entropic contribution in the exponent of D. 
accounts for the elimination of degrees of freedom associated with the replacement of the original Hilbert space Jif 
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by the space ^ that involves a smaller set of observables, those generated by the operators {M} and their products. 
The non-negativity of the operator F is ensured by that of the matrix F. More concrete interpretations of F_ will be 
given in Sec. lVIIll bv Eas. (ll39l) or (1 14411 . 

While (Ma)map = IS evident, the surmise (11081) is proved in Appendix B. It is first shown there that the Kubo 
correlations of the operators {M} in the state F are given by 

^ • ( 112 ) 

\ / map 

The ordinary correlation matrix of the operators {M} is then derived from (111211 and shown to coincide with the 
matrix B defined by (IMll : 

«M<. - = *“-> ■ <‘13) 

Thus, the operator J2 can be regarded as a substitute in to the exact state D oc exp(—/3K) that satisfies the 
anticipated identities (11071) and (11081) : We can interpret the matrix elements of B as exact correlations of the operators 
{JM} in the state D. 

At first order in the sources, the variational treatment amounted to replace the operators {M} by their expectation 
values Here, at second order, we reproduce the variational approximation B for correlations by keeping 

contributions of lowest order in the deviation {M~ both as regards the reduction (110911 of the algebra of {M} 

into that of {M} and as regards the expansion (IllOl) of the free energy operator. 


B. Heisenberg dynamics of the mapped Lie algebra 


Let us extend the above results to the time-dependent correlation functions. We restrict here to the case where 
ff = K, as in Sec. lVl Al As already seen, the solution of the approximate backward Heisenberg equation is then 
generated by the kernel jCF according to Eas. dMll for the coordinates a ^ 0 [and to Eg. (1701) for a = 0]. This 
time-dependence keeps the Heisenberg operators in the Lie algebra {M} since they are given by 


or equivalently by 








dt 


(CF)f 


(114) 


(115) 


The equations (|115|1 constitute a linear set. We are thus led to define, in the mapped space the time-dependent 
operators Ma(^^0 by the corresponding equation 


= -(CF)i^ (Mg(t',t) - Rf ), (116) 

with the boundary condition Ma(^!0 = Ma- Moreover, from the definition (|111|1 of the operator F and from the 
mapped algebra (11091) . it follows that 


[M„,£] = *?i(CF)i'(M^-i?f). 


(117) 


Hence the equations of motion of the set {M^(^^0} defined by (11161) read alternatively 


d&,t) 

dt 


1 

i h 




(118) 


One recognizes the structure of an exact backward Heisenberg equation in the space where the free-energy 
operator F plays the role of a Hamiltonian. As in the static case, the mapping of the operators {M} onto {J\/[} replaces 
approximate properties in Jif by exact ones in here the approximate dynamical equation (11151) in the space Jff 
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by (|118l) where the commutator is restored in the space Accordingly, the mapped Heisenberg operators 
are given by 


(119) 

a mere unitary transformation in the space whereas the transformation (11141) of the set {M} in the space is 
not unitary. 

The linearity of Eas. dllbL or equivalently the quadratic nature of F, exhibits harmonic-oscillator dynamics in the 
mapped space. This will be made more precise in Secs. IVlll PI and IVlll El In Eg. (11181) the effective ’’Hamiltonian” F 
in the space should not be confused with an approximation for the Hamiltonian H = K oi the original problem. 
Its expressions (lllip and ([56|) are related to the free-energy function f{R} and its deviations, rather than to the 
original Hamiltonian H. The contribution to F of the entropic term —through F = K — /3“^S, is essential. This 
contribution is the only one left if K belongs to the Lie algebra, and it remains finite at zero temperature. 


C. Unified formulation of the variational expressions 


Our mapping has provided a formalism ('Sec. lVH Al) in which both optimized expectation values and correlations 
in the state D are generated, in the mapped space as traces over the effective density operator F oc exp(—/3F). 
This operator depends on the temperature both explicitly and through the kernel F of the operator F. 

An arbitrary observable Qj acting in the original Hilbert space is now represented by its mapped image given, 
according to (ESI), by 


Q^. = Mo + (M„ - i?® Mo) ^ 0) • (120) 

uRot 

Then, the optimized expectation value of a single operator Qj is 

and that of the product Qj Qfe is 

(Q. Ofe)app = (122) 


both being directly obtained from the effective state D_ in the space . 

In the case H = K considered in Secs. IVI Al and IVH B[ the time-dependence of the Heisenberg operator associated 
with Q^. is governed by the backward Heisenberg equation (I118|) , which yields 


Qj{t) = . 

Hence, the variational two-time correlation functions take in the mapped space the simple form 

C,,{t',n = {TQf{t'-t")Q^) -(Q) (Q) . 

—J —^ map —J map — 3 map 


(123) 


(124) 


The same operator F [Eg. (11111) ] occurs both in the density operator F and as an effective Hamiltonian in (11231) . It 
depends on the following ingredients: R^a'^ defined self-consistently by Eq. EJ), Ca/ 3 {F(°)} given by T'fp R^\ F given 
by (ISSl) and D by (|llll) . 


VIII. THE EIGENMODES AND THEIR INTERPRETATION 

In this section, we consider only initial states at equilibrium, in which case the dynamics is governed by the 
Hamiltonian H = K. The variational expressions of correlation functions obtained in Sec. |V] then involve functions of 
the sole matrices C and F. Practical evaluations rely on their diagonalization. The eigenvalues and eigenvectors thus 
obtained will enable us to interpret the above results. 





























25 


A. Diagonalization of the evolution kernel and the correlation matrix 

The product iCF of the commutation matrix C and the stability matrix F is the kernel which governs the evolution, 
as exhibited by the dynamical equations (IMl) of Sec. IVI Al It also occurs in the expression (|55)) of the correlation 
matrix B. The diagonalization of the matrix iCF, and the study of its properties, appear therefore appropriate. 

The Lie algebra is globally hermitian, namely, if the operators of the basis are not individually hermitian, they 
come in conjugate pairs (e.g., and in the fermionic example). We denote as a the index of the operator 

Mq = MJ,, which may or may not differ from Mq,. The change a i—>■ a of all the indices in the quantities Ra, Q“, Cap 
or F“^ transforms them into their complex conjugates. Hence, iC, which is antisymmetric, and F, which is symmetric, 
are equivalent to hermitian matrices. Moreover, F is non negative. (Only the case of a strictly positive matrix F is 
examined here; vanishing eigenvalues are considered in the end of Sec. lIX Cl l Altogether, the matrix iCF is equivalent 
to the antisymmetric hermitian matrix iF^^^CF^/^, so that its right and left eigenvectors, respectively denoted as ip 
and (j), constitute a complete biorthonormal basis in the space a, while its eigenvalues are real and either vanish or 
come in opposite pairs. 

Using the above properties, one can classify into three subsets the eigenvectors of iCF, denoted with indices n, —n 
and p, respectively: 

(i) To the positive eigenvalues > 0 are associated the right and left eigenvectors ?/>" and 4>n defined by (a, /3 ^ 0) 

= (125) 

C * F^'^ = = F^'^. (126) 

(ii) Taking the complex conjugate of (11251) provides the right and left eigenvectors '0”"' and associated with 
the negative eigenvalue — 


c” = mr , Cn = (o* ■ ( 127 ) 

(iii) The eigenvectors and (j)p associated with vanishing eigenvalues are given by 

iCa^¥'^^rp = 0, cj)‘^Cap = 0, (128) 

Va = m* , 4 = (4)* = • (129) 

The biorthonormality of the sets {0} and {0} is equivalent to the orthonormalization relations for the right eigen¬ 
vectors {0} expressed by 


€ rp = m* rp = < 5 ””', 0 ^ 0 ^' = m* 0 ^' = , 

02 F“/5 rp = c F“^ rp = mr f“^ 0^ = 0. (lao) 

Likewise, the closure property is equivalent to 

rjT + mr rp]+Y.^irp- ( 131 ) 

n p 


The matrix iCF is thus diagonalized as 

i{C¥)i = ^U„[0S0^ - (0A)* (0F)*] , (132) 

n 

a sum involving only the modes n associated with > 0. The matrices F and C can then be expressed in the form 

F“^ = E[('('-)* + C (<(-1)1 + E € (133) 

n p 

iCap = Y^ u„[02 (0f )* - m* rp] , (134) 

n 

which for F includes also terms associated with the vanishing eigenvalues of iCF. 

The above diagonalization appears as a generalization, for an arbitrary Lie group and at non-zero temperature, 
of the standard diagonalization of the RPA matrix for fermion systems [H, Is^. The reality of the eigenvalues 
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which follows, as shown above, from the positivity of the matrix F is well known in that case [35l - l37ll. The present 
diagonalization is also similar to the diagonalization of a quadratic Hamiltonian of boson operators [34l| . 

Inserting the expansions Ijl33ll and (113411 of F and C into the expression (IM)) of the correlation matrix B yields 


^q:/3 


= E«r 




mr 




— exp(—/3 0 

li ^ 1 


(135) 


exp(/3 /iH„) — 1 




This expression involves only the eigenvalues and right eigenvectors of iCF, normalized according to (11301) . 
In the zero-temperature limit, only the first term of (jl35l) survives, and B reduces to 


(136) 


where the parameters and are found from the limit /3 —> oo of iCF. The quasi-scalars do not contribute. 

In the high-temperature limit or in the classical limit, B = (/3F)~ is given in diagonalized form by (113111 . The same 
holds for the Kubo correlations. 


B. Diagonalization of the effective free-energy operator 


The above diagonalization of the matrix F will help us to give an interpretation of the effective Hamiltonian £ 
defined in the mapped space ^ [Ea. dlllll of Sec. lVH A] . Let us introduce, besides the unit operator Mo = ij the 
following new basis for the mapped algebra {M} (a 0): 




A — 

5 L^-n —n ’ 


y/Kn. 

Y, = C (Me - = (Yjt. 


(137) 


The pair of operators A„, A_^ is associated with each mode r2„, and the single operator Yp with each vanishing 
eigenvalue of fCF. Conversely, the original operators of the mapped Lie algebra are decomposed over the modes n 
and p according to 


Mo - = E [C A„ + (V'S)* At ] + E V'S Yp/. 

n p 


(138) 


The eigenvectors ip appear as the amplitudes of Ma on the different modes. 

The free-energy operator E takes in the new basis the diagonal form 

F = E i (At A„ -f A„ At ) + E ^ Y^ . 

n P ^ 


(139) 


The commutation relations of the operators A„, AJ^ and Yp follow from those ([Moj Ma] = ihCa/s) of the opera¬ 
tors {Ml; from the diagonal form (I134p of the commutation matrix Ca /3 and from the biorthogonality (11301) of the 
amplitudes ip and (p. The transformation (11371) thus implies 


[A„, Ap/] — Snn' , (140) 

[A„, A„,] = [Ap, Yp] = [Yp, Yp,] = 0. 


C. A bosonic and scalar algebra 

The mapped algebra {M} is spanned, according to (I138L by two sets of operators. On the one hand, the commuta¬ 
tion relations (I140p express that this algebra has a symplectic sector n: the set {A, A^} is simply an algebra of bosonic 
operators, with single-boson states labelled by the index n. On the other hand, in the sector p, the set {Y} commute 
with all other operators of the mapped algebra, and hence each Yp can be regarded as a scalar random variable. 
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The form (|139|) of the free-energy operator F in the space is suggestive. It can be identified with a Hamiltonian 
of non-interacting bosons, with single-particle states n having the energy plus a classical quadratic part with a 

coefficient 1/2/3 for each variable Y^. 

The evaluation of expectation values (•■•)map involves traces over the space ^ with the weight F oc exp(—/3F). As 
regards the bosonic part of (11391) , ^ contains a Fock space and the trace is as usual a summation over the occupation 
numbers of each single-boson state. As regards the scalar part, the trace is meant as the integration ]/[p f dYp. Thus, 
the expectation values over F oc exp(—/3F) of single operators A^j, or Yp vanish, while we find for pairs 

(a). An),nap = ^ . (1^1) 

(Yp*)„.p = (!yg).„ap = i; 

all other expectation values of pairs vanish. Each mode n yields a Bose factor associated in canonical equilibrium 
with the energy no chemical potential occurs for the distribution F oc exp(—/3F). For the scalar variables, the 

variance 1 arising from the Gaussian weight exp(—Yp^/2) within F provides an expectation value 1/2/3 for the term 
Yp^/2/3 of F, in agreement with the equipartition theorem of classical statistical mechanics. 

We are considering in this section the dynamics (11191) for H = K, and we have seen that F then plays the role of 
a Hamiltonian in the space p^. The diagonal form (|139l) of F indicates that the mapped operators A„ oscillate as 

(142) 

while the scalars Yp remain constant. Thus, the operators Ajj,A)j behave as bosonic annihilation and creation 
operators in all respects: commutation relations (|140l) . average occupation in canonical equilibrium (11411) and dynamics 

dmi). 


D. Interpretation in terms of oscillators 

The mapped bosonic Fock space can equivalently be regarded as a space of oscillators, each single bosonic state n 
corresponding to an oscillator mode with frequency fin- The operators A„, A)j are thus replaced by the position and 
momentum operators 

X„ = ^A72(A„ + At), P„ = x/ft72(A„-A|,)/z, (143) 

satisfying the canonical commutation relations [X„, = ihSnn'- The effective Hamiltonian F then takes the form 

f = Y.^n (P„^ + X„^) + ^ E • (144) 

It now describes uncoupled harmonic oscillators, plus a quadratic energy associated with the classical variables Yp. 
In this alternative interpretation, we have 

(x;)„.p = {PS)„.p = 5coth^^, (145) 

and the dynamics describes harmonic oscillators with frequency 


E. Quasi-bosons and quasi-scalars 

From now on, we return to the original space The correspondence between the algebras {M} and {M.} leads 
us to perform on the original Lie algebra the same transformation as (|137l) - (ll38p . which introduces for a given 
a new basis A„, A),, Yp (besides Mq) in the algebra {M}. 

The commutators between these operators, issued from the Lie structure (ED, are not simple. However, the 
expectation values of these commutators, evaluated as traces over have the same structure as the commutators 
(I140p of the corresponding mapped operators A„, A)), Yp, that is, 

([An, A|j,])^pp = 6 nn ' , 

([A„, A„/])jjpp = ([An, Yp])g^pp = ([Yp, Yp/])g^pp = 0. 


( 146 ) 
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The operators A„, can thus be termed ’’quasi-boson” annihilation and creation operators, while the operators Yp 
can be termed ” quasi-scalars”. 

The matrix ihCap was defined as the expectation value ([Mq,, M/ 3 ])^pp in the original basis of the Lie algebra. 
Going to the new basis amounts to diagonalize C according to (|134|1 . In this new basis, iAC takes a form involving 

only 2x2 diagonal blocks ^ q ^ ^^.d zeros, in agreement with the above relations. 

The variational expectation values in the state of the operators A„, Aj^ and Yp vanish, 

(An)app = (ADapp=(Yp)app=0, (147) 

while (Mo) 3 ^pp = 1. The approximate expectation values of pairs of operators, 

(A,\ A4,pp = = (A„ Al),pp - 1, (148) 

(YApp = (MXp = 1> 

are the same as the equations (11411) which were exact in the space again all other expectation values of pairs 
vanish. We recognize, as in the mapped space, the Bose factor of a mode n for the quasi-boson operators A„, A^, and 
the unit variance for the Gaussian quasi-scalar variables Yp. 

The change of basis {M} i—^ {A, Af,Y} is issued from the diagonalization of Sec. lVIII~S1 This change leads to 
re-express the correlation matrix B as 


= (M„ M^)^pp - Rf = ^C (A„ At ),pp (V>f)’ 

n 

+ ^ (^^)* (At A„)^pp rp+^Y. V’S(Y/),ppV’^ 


(149) 


The Bose-like factors that appeared in the expression (11351) of B now come out as expectation values (11481) of pairs 
of quasi-boson operators A„ and At,, while the last term of (11351) is consistent with the variance equal to 1 of the 

quasi-scalar operators Yp. The coefficients arise from the change of basis (11381) that expresses the operators — Ra^ 
on the basis {A, At,Y}. 

More generally, to write the optimized correlation of arbitrary observables Qj and Qk, one should express their 
images on the basis A, At, Y according to 


Q, = Mo + A„ + {Q]^r A],] + ^ S 


iP Y 
■js ' P ’ 


(150) 


where the new, bosonic and scalar, coordinates are (a 0) 

(we recall that = dqj{R^°^/dRa^). One thus finds for the correlations 

(QjQk)a.vD ~ (*5i)ann(Qfe)app = 


= E 


app \ J / app \ app 

p hD,. 


coth 


ReC"b(«b) +ninC"i.(«b) 




P 

ks ' 


The antisymmetric part of the correlations, 


1 


p:{QjQk - QfeQ,)app = (Qfcb)* , 


(151) 


(152) 


(153) 
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agrees with [Qj, Q^] =ih Cq /3 Qf An alternative interpretation of the symmetric part is obtained by 

replacing, as in Sec. lVIIlPI quasi-bosons by quasi-oscillators with variables X„, P„. Writing the image of Qj in the 
basis {X, P, Y} then produces the factors (S/2) coth(/3 Sn„/2) which are the expectation values = (Pn)app' 

Since H = K here, the approximate Heisenberg operators in the space follow the same evolution as (11421) in the 
mapped space that is, 


t) = t) = Yp . 

Hence, the two-time correlation function Cjk(t',t") is given, when the initial state is at equilibrium, by 

g - i 


Cjkit'X) = 




(Q/b) 


it) ^ ^ — ^ h fir. 


(Qfcb)* 


(154) 


(155) 


g/3Rn„ _ 1 


Q 


kh 


J2a%ai. 


which is the explicit form for the general equation (11241) . It exhibits the coefficients of the images Qj and Qk on the 
quasi-boson and quasi-scalar basis, the Bose factors and the boson dynamics. 


F. Linear response and excitation energies 

The linear response (llOOl) is directly found from (IIOII) and (11501) : it does not involve the quasi-scalar contributions 
nor the Bose factor. Its dissipative part, defined through a Fourier transform with respect to t' — t”, comes out as 

J E <5(^ - (Qfeb)* - {Q]S S{u: + H„) Ql^] . (156) 

n 

Thus, the present approximation yields the H„’s as resonance frequencies. Consistency properties such as the Kramers- 
Kronig dispersion relations and the Kubo fluctuation-dissipation relations are satisfied. 

At the zero-temperature limit, the exact expression of Xjki^) given, in terms of the ground state |0) of K and of 
the excited states |exc) with excitation energies Aexc, as 

Xjki.^) = 7r^[(0|(5j|exc)(exc|Qfe|0) 5{hw - (157) 

exc 

-(OlQfelexc) (excIQjlO) 5{huj -I- F^exc)] ■ 

Comparison with (11561) shows that the positive eigenvalues of i?iCF can be identified as approximations for the 
excitation energies En of some set of states, labelled as \n). The amplitude Q/j, of the image Qj over the quasi-boson 
annihilation operator A„ appears as an approximation for the matrix element {Q\Qj\n). In particular, the quasi-boson 
operators A„ satisfy approximately 


(0|A„|n') = 5nn' , 

as if the state \n) were obtained by creating a quasi-boson. 


(158) 


IX. SMALL DEVIATIONS 

We consider in this section static and dynamic changes brought in by a shift of the initial conditions. 


A. Static deviations 


The change in the equilibrium properties arising from a variation 5K of the operator K that defines the state 
I)(xexp{—PK) is accounted for in the variational treatment by the shift of parametrized by the set 


= Tr5P(°) M« = Tr[p(°) -h (5P(°)](M„ - i?®). 


(159) 
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Denoting by 5K. the image of 5K^ the first-order change in the self-consistent equations o provides 

= /3 (5/C“ -F /3 5Rf , (160) 

that is, using the definition (1551) of the matrix F, 

= (a,/3 7^0). (161) 

The resulting change AF of the free energy F ~ min/{i?} = min[fc{i?} — T 5'{i?}], expanded up to second order in 
SK, is found as 

AF « i SR^°^ SR^p^ + SIC°‘ SR^°^ (162) 

= (5fc - i 6R‘^'> F“^ Sr‘'^^ = 5k- i(5/C“ (F-^)„^ (5/C^ . 

In (|160|) - (|162p . the matrices K, S, F, the image 51C and the symbol 8k of 8K are taken at In particular, the 

expression (153 of the heat capacity is recovered by taking 5K = K5/3//3. 

This shift takes a suggestive form if F“^ is replaced by its diagonalized form (11331) . Using (11591) and (I147L and 
denoting as 


= Tr -t 5V^°^] A„ = Tr 5V‘'°^ A„ , (163) 

5R^^* = Tr 5T>(°^ Al, 5R^°J = Tr Yp 

the variations of the symbols of the quasi-boson and quasi-scalar operators A„, A)j and Yp, one finds from the variation 
of f{R} the form 


_ i SR^°^ F“^ SR^ (164) 

= 5k-Y,hn^ ^ ^ {5R^^f . 

n p ^ 


( 0 ) ^ 

Each oscillator mode provides a contribution hQ.n weighted by the amplitude , while each quasi-scalar mode 

brings in the equipartition contribution (2/3)”^. 

B. Dynamic deviations 

Let us turn to the change in the time-dependent expectation values induced by a shift around 

the equilibrium state The resulting small deviations dRa\t) around Ra\t) are governed by the set of equations 

= (L + C H)f SRf {t), (165) 

obtained by varying Ra\t) in Ea. (l53 : the initial conditions 5Ra\ti) are given by (11591) . [For the fermionic single¬ 
particle Lie-algebra, p65l) is the dynamical RPA equation issued from the TDHF equation.] The kernel L -|- CH 
governing the forward equation (|165l) for the small deviations turns out to be the dual of the kernel governing the 
backward equation dfiZl) for the approximate Heisenberg observables The shift 5{Qj)^, variationally given 

by S{Qj)^ = Q'^{R‘^°\t)}5Ra\t), is therefore equal to 

S{Qj),c^Qf^{t,nSR^°Ht") 

for any intermediate time t”. This independence on t" agrees with the expression dZOD of {Qj)^. 


( 166 ) 
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C. Dynamic and static stability 


We now specialize to the case H = K, for which an initial shift {SR^^^} generates a deviation {SR^^\t)} of 
around the fixed value {R^^^}. The linearized equations (11651) then involve the constant kernel L + CH = CF, and 
can be solved as 




, (t-ti)CF 


SR 


(0) 


(167) 


Here, as in the backward equations (I114I) . (I115I) . the dynamics is governed by the product CF, so that the result (11671) 
can alternatively be written as 


5R^^\t) 


, (i-ti) CF 




(168) 


Within the variational treatment we recover for small deviations the equivalence between Schrodinger and Heisenberg 
pictures, already exhibited in Eq. dZQl). 

If the matrix F is positive, one can change the basis so as to diagonalize iC F according to (jl32l) and use the new 
variables SR^j^ = TrSV^^^A„ and (5i?ps^ = Tr SV^^Wp, which yields 


= ( 169 ) 

SR(^J(t) = SR(°J . 


These evolutions merely reflect the motion (jl54l) of the Heisenberg operators {A^, Y^}. Again the modes are decou¬ 
pled, quasi-bosons (or oscillators) undergo pure oscillations and quasi-scalars are static. 

The above sinusoidal form of the dynamics entails Lyapunov stability. This property is defined, according to [s^, 
as follows: ’’The equilibrium point xg is said to be Lyapunov stable if given any neighborhood U of xo, there is a sub¬ 
neighborhood V of Xg such that if x lies in V then its orbit remains in U forever.” In other words, the trajectories tend 
uniformly to as their initial point {R^^^ti)} tends to Hence, the Lyapunov stability of the linearized 

motion is ensured if all eigenvalues of F are positive, that is, if the approximate free energy associated with is 
a local minimum of f{R}- Such a property is well known for fermions, both at zero [s^ and at finite temperature 
in which case the minimization of the Hartree-Fock (free) energy entails the reality of the RPA modes. We 
find here, in a general variational context for any trial Lie group, that the static stability of the approximate ’’state” 
implies the dynamic stability of motions SR^^^t) around it. 

The matrix F may have vanishing eigenvalues. This occurs, for instance, if a continuous invariance is broken by 
the approximation for the equilibrium state; in this case f{R} is minimum for a continuous set of solutions 
Since F is then not invertible, iF^/^CF^/^ is no longer defined and it is not ascertained that the matrix *CF 
is diagonalizable (some right eigenvectors may be missing). If iCF is diagonalizable, all its vanishing eigenvalues 
yield a constant contribution to the set SR^^^t) so that Lyapunov stability is still ensured. However, if it is not, 
contributions behaving as powers of t come out, so that the dynamics is unstable (for a more detailed discussion, see 
Appendix B of [1^). Such a behavior may be associated with Goldstone modes. For instance, if the original problem 
is translationally invariant and if this invariance is broken by a localized solution a small shift in the initial 

conditions produces an instability, characterized by a boost at some constant velocity. 


X. CLASSICAL STRUCTURE OF LARGE AMPLITUDE MOTION 

We have written in Sec. lIVBl several equivalent dynamical equations for the time-dependent expectation values 
{R^^^{t)} of the operators {M}. We now analyse the structure of these equations, which will turn out to have a 
classical form for any Lie group. 


A. Poisson structure 

Let us consider {R} as a set of classical dynamical variables, and let us show that the tensor <Cap{R} = Lj^A^ 
can be regarded as the generator for this set of a Poisson structure issued from the Lie algebra (ED. We first recall 
the definition of Poisson structures (40l l4lj| . Consider some set of dynamical variables {x} parametrizing points on 
a manifold, and functions /, g, h, ... of these variables. A Poisson structure is a mapping .f/, ^ from a pair of 

functions f,gtoa third function k, which obeys the following rules: 
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(i) bilinearity; 

(ii) antisymmetry: ff,g^ = 

(iii) Jacobi identity: H/, gj, kj + Hfc, f}, g} + k}, kfj = 0; 

(iv) Leibniz derivation rule: ifg, kj = ffg, kj + ff, kjg . 

We consider here functions g{R} of the variables {R}, and define a Poisson structure through 

RpJ = C^p{R} = TlpR^ (170) 

and 

(m) 

One can readily check that the above rules are satisfied, owing in particular to the properties of the structure 
constants namely, antisymmetry and Jacobi identity. The Poisson structure (I170D - (I171I) is thus generated by the 
Lie algebra of the set { M }. 

The equations of motion (1621) for can now be rewritten as 

= h{R(°Ht)n. (172) 

These quantum variational equations are therefore identified with classical equations involving the brackets (11711) and 
governed by a classical Hamiltonian h{R^^^ (t)}, the symbol of the quantum Hamiltonian H. 

The relation (ICT) expresses the time-dependent expectation value of the observable {Qj)t as a function of {i?^°^(t)}; 
together with (11721) . it implies that {Qj)t evolves according again to the classical dynamical equation 

= iq,{R^°Ht)}, ft{i?(°)(t)}:& , (173) 

which involves the symbols of Qj and H. The variational approach, together with the introduction of symbols, thus 
generate approximate dynamics of expectation values that have a classical structure, generated by the Lie-Poisson 
bracket p71D and by a Hamiltonian. 

The symbol h{R^^\t)} is obviously a constant of the motion. Moreover, the von Neumann entropy S{R^°\t)} 
defined by ((^ is also a constant of the motion. Indeed, using (l26t then (l35l) . one finds 

iRf\t), = 0, (174) 

(J yi) 

which implies that S{R^^\t)} remains constant during the evolution (11721) of {R^^^{t)}. 

Lie-Poisson structures for dynamical equations issued from a variation al p rinciple have been recognized in cases 
such as the Vlasov equation [4^, time-dependent Hartree-Fock equations [^, time-dependent Hartree equations for 
bosons [l2| and for (j)^ field theory [l^, [l^ . 

Proposals of non-linear extensions of quantum mechanics have suggested a formulation in terms of a Poisson 
structure . Here it is the restriction of the algebra of observables to the trial Lie algebra which produces a Poisson 
structure within standard quantum mechanics. 

B. Canonical variables 

We have seen in Secs. lVHI HI and fVHI PI that the diagonalization of CF generates in the mapped space a linear 
transformation (jl37l) . (jl43l) of the operators {M}, which produces pairs of canonically conjugate operators X„, P„ and 
scalars Y^. This corresponds in the original Hilbert space to a construction of quasi-oscillator operators X„, P„ and 
quasi-scalars Yp (Sec. lVIlf^ . Accordingly, the expectation values of their small deviations, defined by 

JA„(t) = (t)] = ^/^Tr JP(°)(<) (A„ + At), 

JP„(t) = (i) - = vW2Tr JpW(t) (A„ - At )/*, 

evolve according to (11691) as conjugate variables of classical harmonic operators with frequency f2„, while the quantities 
SYp = SR^pJ = TTSV^°\t)Yp remain constant since Fp, = 0 for any h. A classical symplectic structure thus 






























33 


appears for linearized motions, with ordinary Poisson brackets {6Xn,SPn>} = Snn', {SXn,SXni} = {6Pn,SPn>} = 0 
and Hamiltonian (1/2) + X^^) . 

The above canonical classical structure pertains to the dynamics of small amplitude motions of + 

around equilibrium We have seen however ISec. lX All that a more elaborate Poisson structure occurs 

for large amplitude motions. In this case, the bracket ^Ra,R/ 3 ^ = Cq/jII?} depends on the dynamical variables, 
whereas for linearized motions we had simply to diagonalize the constant matrix according to (11341) . 

Moreover, the motion is no longer harmonic. 

Nevertheless, one can re-express the classical dynamics of Sec. IX Al in terms of canonical variables by means of a 
non-linear change of variables in the space {R}. Indeed, a ’’splitting theorem” states than an arbitrary Poisson 
structure can locally be split into symplectic components and invariant components: There exist independent (non¬ 
linear) functions of the coordinates {i?} such that their brackets (11711) reduce to 

— djin' j ^n' '^n' (1^^) 

Vp^ ~ Vp^ ~ ^Vpt Vp'^ ~ 0 ■ 

In general such a transformation does not exist globally but only locally in some neighborhood around each point of 
the trajectory of the point {i?^°^(f)} in the space {i?}. The dynamical variables and 7r„ are canonically conjugate 
in the elementary sense (their brackets reduce to ordinary Poisson brackets). These variables are dynamically 
coupled and their motion (11731) is governed by Hamilton’s equations, while the variables Pp (the Casimir invariants) 
are structurally conserved in the flow. 

The construction of the set {^, tt, rj} is not simple and does not result from the mere diagonalization of CapiR} at 
each point, contrary to the linearized dynamics. In fact, the reduction of the Lie-Poisson structure to the form (11751) 
has been achieved only in special cases; relevant to the present work are the Vlasov equation for which a rigorous 
proof has been given [i^, the time-dependent Hartree-Fock theory at zero (T^. l43 - l4^ and finite temperatures. 


C. Stability of equilibrium and of non-linearized motions 

In Sec. lIXCl we have shown that small amplitude motions around thermodynamic equilibrium {7?^°^} governed by 
the linearized equations for {R^^^{t)} (with H = K) are Lyapunov stable if the trial free-energy function f{R} is 
minimum at {7?^'^^}. Let us show that this stability property also holds for motions around equilibrium governed by 
the non-linearized equations of motion (1621) . 

To this aim, we rely on the classical form (11721) of this equation written in the Poisson formalism. We note moreover 
that, owing to the property (11741) of the entropy function S{R^°\t)}, the addition of —/3“^S'{77(°)(t)} to the classical 
Hamiltonian h{R^^\t)} = k{R^^\t)} does not affect the dynamics, so that the free-energy function also generates the 
large amplitude trajectories of {R^^\t)}: 

^ C„47?(°)(7)} . (176) 

(t) 

The time-dependence of f{R^°\t)} is given, according to (I173|) . by df{R^°\t)}/dt = ff{R^°\t)}, f{R^°\t)}J = 0. 
Hence f{R^^\t)} remains constant along any trajectory. On the other hand, if {7?^°^} is a stable static equilibrium 
point, /{7?} has an isolated minimum equal to /|77^°H = F- These two conditions ensure Lyapunov stability according 
to the Lagrange-Dirichlet theorem (see for instance [^): the trajectory {R^^\t)} uniformly tends to the point 
if the initial value f{R^^\ti)} tends to F. 

Thus, for non-linearized as well as for linearized motions of {7?*^°)(t)}, the present variational approximations 
preserve the following property of the exact dynamics near an equilibrium state: the static stability, f{R} minimum 
at {7?} = {7?(°^}, entails the Lyapunov stability of the dynamics in some neighborhood. 


XI. RESUME OE OUTCOMES AND DIRECTIONS FOR USE 

We have dwelt at length on the derivation of the results of a variational approach based on the principles of Sec.HIl 
and on the restriction of trial spaces to Lie groups. We summarize below some of the formal outcomes thus obtained; 
they are sufficient for practical applications to specific many-body problems. 
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A. General features 

The initial step consists in selecting, among the whole set of operators in Hilbert space, a Lie algebra spanned by a 
set of operators {M} labeled by the index a. (We include in this Lie algebra the unit operator I, denoted as Mq.) We 
gave above as seminal example a system of fermions for which the set {M} encompasses the single-particle operators 
aj^ciiy, the index a denoting the pair The approach can be applied to many other systems. For fermions with 

pairing, the Lie algebra {M} includes the operators and alaj^ (with fi > u); here the index a denotes not only the 
pair but also distinguishes between the operators aj^ajy, and alaj^. For bosons with possible condensation, 

one can take as basis {M} for the Lie algebra the operators I, a^, aj^, a^ai, {fj, > v) and (/r > i/); coherent 

states arise from the inclusion of the operators and a|^. If such a bosonic system is translationally invariant, the 

algebra can be reduced to I, gq, aj, aj^a^, akG-k, (where k denotes the momentum of single-particle states). 

The formalism also applies to other sub-algebras, or to other systems such as mixtures of fermions and bosons, spin 
systems or quantum fields as in [l^ , the only restriction being the Lie algebraic structure of the set {M}. In all such 
cases the algebra {M} is characterized by the structure constants entering the commutation relations 

(177) 

As we wish to derive static or dynamic properties of some observables Qj of interest, we have relied on the variational 
evaluation of a generating functional 1 Secs. UTAI and HTCl) . so as to deal simultaneously and consistently with all such 
properties. This has entailed the introduction of two sets of trial objects, namely, trial ’’generating operators” A that 
depend on the observables Qj and their associated (time-dependent) sources, and trial density operators V. Both A 
and V are elements of the Lie group generated by the chosen Lie algebra {M}, that is, exponentials of elements of this 
algebra. We focused on expectation values and pair correlations of the observables Qj; they are found by expansion 
of the generating functional in powers of the sources. We summarize below the results thus obtained. 

We parametrize the trial operators T), which behave as non-normalized density operators, by their normalization 
Z = TiV and by the numbers Ra = TrlVIaP associated with the operators Mq (a 0), where V denotes the 
normalized operator Z~^'D. The exponential form of the operator "D is suited to investigate systems at non-zero 
temperature; ground state problems are dealt with by taking a zero-temperature limit. We may thus approximately 
answer questions about an equilibrium (unnormalized) state D = exp(—/SAT); for a system in canonical equilibrium, 
K is the Hamiltonian. Other choices of K allow us to deal with non-equilibrium problems, D = exp(—/SAT) being 
then the initial state. 

An essential tool consists in the representation of the observables Qj, of the operator K entering D = exp(—/JAT) 
and of the Hamiltonian H by their symbols qj{R}, k{R} and h{R} within the Lie group (Sec. lHI SI) . Defined for any 
operator Q and an arbitrary element T) of the Lie group as 

q{R} = TYQf>, (178) 

a symbol is a function of the scalar variables Ra {ex ^ 0) that parametrize R. The practical implementation of 
the approach is based on the possibility of evaluating explicitly such symbols ( 7 {A}. (In the fermionic and bosonic 
examples, this is feasible owing to Wick’s theorem.) One also needs to express in terms of the variables Ra the entropy 
function Aji?} = —Tv'Dhx'D. 


B. Static quantities 

The generalized free energy F = —InTr exp( —/SAT) and the thermodynamic quantities issued from it are 
obtained fSec. lIV Al) by looking for the minimum of the trial free-energy function f{R} = k{R} — T ^ji?} (where 
k{R} is the symbol of AT); the values of Ra = Ra'^ at this minimum are given by = 0. The resulting, consistent, 

approximations are then, for the free energy: 

r _ „ 3fcfR(0)} . . 1 

A =-/3“4nTre“'’'^ ~=-/3“4nTreL ^ H, (179) 


for the entropy: 


S-S{R^°'>} = 


9/{A(0)} 


(180) 


dT 
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for the energy (when K is the Hamiltonian): 

(K) ~ , 

and for the expectation values of the observables Qj in the state D oc exp {—jSK) iSec. lIVB|) : 

TrQf 


{Qj) = 




Tre-/3^ 

Thermodynamic coefficients, such as the specific heat [Ea. (l57)) ]. are found from the second derivatives 


(181) 


(182) 


(183) 


which appear as a matrix in the space of indices a. 

While the above expressions appear as mere extensions to arbitrary Lie groups of standard mean-field results, the 
approach has also provided variational expressions for the correlations in the state D oc exp {—PK) (Secs. |VD] and 
IV E|) . For correlations between the elements of the Lie algebra, we have obtained {a,p ^ 0 ) 

(M„M,)-(MJ(M,, = (184) 


where the matrix C is defined by . More generally, the correlations between two observables Qj and 

Qk have been found as 


{QjQk) - (Qj) (Qk) 


dqk{R^°^ 

dR‘'°^ 


^a/3 


dqk{R^°^ 


(185) 


in terms of the matrix B and of the symbols of the operators Qj and Qk- Fluctuations AQj result for j = k. In the 
classical or high-temperature limit, B reduces to (/3F)~^. 


C. Time-dependent quantities 


Here D oc exp {—PK) is the exact initial state at the time ti, and the subsequent evolution is governed by the 
Hamiltonian H. The above results hold as variational approximations at the time p. For later times, the optimization 
of the expectation value {Qj)^ yields (Sec. llVB|) 


iQj)t — 


Tr e ~ t/h qS ^ — i h tjh 




where the time-dependent expectation values Ra\t) of the operators are given by the equations 


dt 


(186) 


(187) 


with the initial conditions RaHu) = Ra^. The effective Liouvillian L is expressed self-consistently in terms of the 
symbol h{R} of the Hamiltonian by 


Lj{i?} = r 


7 


dh{R} 

dRp 


(188) 


(In the presence of non-vanishing structure constants T^^, R^\t) should be replaced by 1 in the term 7 = 0 of (I187|) .l 
We have shown (Sec.|X]) that the dynamical equations (11871) and the resulting ones for {Qj)^ have a classical 
structure, where CapiR} = Rj appears as the Lie-Poisson tensor associated with the classical coordinates {R} 
and where h{R} behaves as a classical Hamiltonian. 

Here again, non standard variational results have been obtained (Sec. IV Cl) for two-time correlation functions. Their 
expression involves an approximation (Sec. IV AT) for the observables in the Heisenberg picture defined by 


Qf qS g-i H(t'-t)/h ^ 


(189) 
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where t' is the usual final running time and t is a reference time at which Q^{t + 0,t) reduces to the Schrddinger 
observable Q®. This approximation ~ ,t)Ma is characterized by the differential equations 

= + (a,/3y^0), (190) 

in terms of the reference time t which runs backward from t' to the initial time ti. The boundary condition is given 

by 


dRS\t') 


(191) 


In Ea. ()190|) . the three matrices L defined by (|188|) . <Cap{R} = ^^apR-y H, which denotes the matrix of second 
derivatives of the symbol h{R\ with respect to the variables i?Q,, are evaluated at the point R = 

The optimized expression for the two-time correlation functions Cjk{t',t") then reads (Secs. |VC] and |VD]) 


L.jk(t ) - Tre-0^ 

^ U) ti) it' > t"), 


{Qj)t'{Qk)t" 


(192) 


which involves both the approximate correlations (11841) at the initial time and the approximate Heisenberg operators 
given by Eqs. (119011 and (I191I1 . 


D. Properties of the results and consequences 

Special cases ('Sec. lVIl) include time-dependent fluctuations, obtained from Cjj{t + 0,<), and linear responses 


X,fc(t',t")= il^mt' - t")[C,kit'R") - Ckjit\t')] 

9it' - t") Qf'^it',t)Co.p{R^°\t)} Ql\t",t ), (193) 

where t is arbitrary in the interval ti < t < t". If the initial state D oc is in equilibrium, with a dynamics 

generated by H equal to K (or to K plus a constant of motion), the approximate expectation values {Qj)t remain 
constant. The equations (11911) are solved as 


Qf'it'A) = Q^{R^°'>} 


gCF(t'-ti) 


Q! 



(194) 


Hence, two-time correlation functions are variationally given in this case by 


Cjkit'.t") 






ihC¥ 


1 — exp(—f hpC¥) 


¥ 


-1 




- CX^ 


(195) 


an expression which, as it should, depends only on the time difference t — t'. Derived within a unified framework, these 
results satisfy consistency properties and conservation laws. For instance, the approximation preserves the relation 
between the static stability of a thermodynamic equilibrium state (/{i?} minimum) and the dynamical stability of 
the motions (jl87l) around it. 

Other special cases are also considered in Sec. ED including classical, high-temperature and zero-temperature limits, 
as well as Kubo correlations. Results for small deviations are presented in Sec. m 

In Secs. IVlJ and IVIH ClTVHI El we have exhibited, for any Lie group, an attractive interpretation for the static 
correlations (11841) and for the two-time functions (I194I) - (I195I) . To this aim, the Lie algebra {M} has been mapped into 
a simpler Lie algebra {M}j tbe operators Ma of which are linear combinations of creation and annihilation bosonic 
operators A„ and (or equivalently of position and momentum operators X„ and of harmonic oscillators), 

and of scalar random variables Yp. Then, the dynamics (11941) is mapped onto a dynamics governed by an effective 
Hamiltonian 

E = ^ 1 (aU„ + A„ At) + E 

n P ^ 

= E^"(^n'+^n') + 4E 

n p 


( 196 ) 
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which describes uncoupled modes. The above correlations (11841) and (11951) take for an arbitrary Lie group the form of 
expectation values of pairs of quasi-boson operators and of scalars in a canonical equilibrium state with Hamiltonian 

F. 

The explicit expression 




E 


Q? 


, — i r2n {t' —t") 


(2?b)* 


I — jShQr. 
p i {i' — 't") 




jbJ g/3fin„ _ 2 


Q 


kh 




(197) 


derived in Sec. IVlllI relies on the diagonalization (11321) of the matrix zCF, which enters the formalism at several 
places. The coefficients and appear as coordinates, in the new basis of the Lie algebra, of the observable 
Qj, their explicit form being given by Eos. (11511) in terms of the eigenvalues ±Hn and eigenvectors of the matrix * CF 
(which satisfy the equations (11301) . (11331) and (11341) 1. Special cases of (11971) are the static correlations written in the 
diagonalized form (11351) and the linear responses, which satisfy the Kramers-Kronig dispersion relations and the Kubo 
fluctuation-dissipation relation. 


XII. CONCLUSION 


We have presented above a variational approach to the determination of various physical quantities pertaining to 
many-body systems. Although our scope has been formal, the generality and flexibility of the treatment appear well 
suited to many specific problems. The results, listed in Sec. EH have been obtained by merging several ingredients: 

(i) The evaluation of a generating functional f Sec. Ill AT) has allowed the simultaneous optimization of different quan¬ 

tities such as the (static or time-dependent) expectation values and correlations of the observables of interest: they 
are obtained by expanding in powers of the sources the functional = lnTrA(ti)Zl. Deriving the correlation 

functions as second-order contributions in the sources of this generating functional provides for them non-trivial ap¬ 
proximations, even for a simple restricted trial space. As an example, for a system of interacting fermions, the use of 
independent-particle trial objects leads to standard mean-field theories for expectation values but to expressions of 
the form (11921) for correlation functions. It is the dependence of the trial objects on the sources which leads to such 
elaborate results. Moreover, approximating all quantities in a unique framework preserves consistency properties. 

(ii) The variational principle for the optimization of the generating functional is built by means of a general method 
fSec. Ill Cl) . The object to be optimized is characterized by some simple equations, regarded as constraints on the 
ingredients, namely, on the initial state D oc exp(—/SAT) and on the time-dependent observables of interest Q^{t',t). 
Lagrange multipliers are then associated with these constraints. 

(iii) For dynamical problems, the Heisenberg observables Q^{t',t) enter the generating functional, together with 
the sources fj{t), through the ’’generating operator” 

A(t) = Te^-^t 

The operator A(t) is characterized by the differential equation (ITT]) expressing , which plays the role of a constraint 
in the variational principle. This Eq. (HU has been derived as a consequence of the backward Heisenberg equation 


dQfjt'A) 

dt 


-[Qfit\t),H] 


(199) 


We recall that, while the standard (forward) Heisenberg equation describes the variation of the Heisenberg observable 
Qj(t',t) as function of the running time t', (11991) is a differential equation with respect to a reference time t which 
runs backward from t' to the initial time A. We have explained in Sec. Ill Bl whv this backward dynamics is the suitable 
one. 

(iv) As regards the initial state D = we have characterized it by the Bloch equation 


dX>{T) 

dr 


FKV{t) =0, 


( 200 ) 


where r runs from 0 to /3. We thus deal with finite temperature. Ground state problems are treated in the limit 
/3 —^ 0; in fact, this procedure turns out to be more convenient than the direct implementation of ground states. 
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(v) The variational equations associated with the above constraints on A{t) and D(t) have been worked out by 
taking a Lie group as trial space for A{t), D{t) and for their associated Lagrange multipliers fSec. lIIIl) . thus replacing 
operators in Hilbert space by functions, their symbols. The approximate states that occur in mean-field theories belong 
to such Lie groups, for instance static and dynamic states in Hartree-Fock approximations for fermions, Hartree-Fock- 
Bogoliubov states for fermions with pairing, or coherent states for bosons. However, the trial operators D used here 
depend on the sources; they are not overall approximations for the exact state D, but serve only to optimize the 
generating functional (and hence all quantities of interest). 

We have recalled in Sec. IXII the outcomes of the above approach. Some well-known approximate results have been 
recovered in special cases, for instance, static and dynamic mean-field or RPA treatments which now appear within 
a general unified framework based on the use of a Lie group as trial space. Moreover, the new variational results 
(I184|) . ()185l) . (11921) and (|195l) have been derived for correlations and fluctuations. In particular, we have shown for 
H = K that quasi-bosons (or quasi-oscillators) come out for any Lie group, and not only in the usual context of 
zero-temperature RPA for fermions. Let us add a few comments. 

In order to determine the generating functional = lnTrA(ti)Il, we have characterized A{ti) and D = e~^^ 

by introducing the functions Aft) and D{t) = e~'^^ determined from the simple differential equations ([Tl]) and ()200|) . 
To build a tractable variational principle, we have been led to introduce the trial time-dependent quantities D{t) and 
Aft) whereas we need only their boundary values I?(/3) and Aft,). Moreover, the number of variables is doubled by the 
introduction of the Lagrange multipliers A{t), associated with the equation for D{t), and V[t), associated with the 
equation for Aft). These quantities are not coupled by the exact equations of motion. However, in a restricted trial 
space, the stationarity conditions entail a coupling between them which allows a better optimization of the generating 
functional 

In spite of their resemblance with density operators and exponentials of observables, the trial objects L){t), A{t), 
A{t), Dfl) are only computational tools for the evaluation of expectation values and correlation functions. In particular 
the trial quantity 'Dft), which looks formally like a density operator in the Schrddinger picture, adapts itself to the 
question asked, namely to the value of '(/'{'f}- If thus depends on the sources (and is not even hermitian). It cannot 
be interpreted as an approximate state of the system, and in fact there exists no approximate density operator in the 
original Hilbert space that would produce both the optimized expectation values (|186l) and correlations (11841) . (11851) . 

(M . 

Since all the results have been derived from the same variational principle, they naturally satisfy some consistency 
properties fulfilled by exact quantities. For instance, if an observable Qj belongs to the Lie algebra and commutes with 
H, exact conservation laws express the constancy in time of the expectation value {Qj)^ and of the fluctuation AQjft). 
These two properties are ensured by the variational approximations ISec. IVI Cl) . In contrast, a naive evaluation of 
fluctuations based on the approximation TiQj^ Dft) ~ TrQ^^ 'D^'^\t) where 'D^'^\t) evolves according to (|5(I)) . would 
provide an unphysical time dependence: The Schrddinger picture for T>A) ft) is variationally suited to the evaluation 
of expectation values, but not of correlation functions or fluctuations, which involve here the approximate Heisenberg 
operators Qf{t, t'). The variational approximation (11951) fulfils another consistency property, as discussed in Sec. IVI Al 
namely that any correlation function Cjkft',t") depends only on the time difference t' — t" when H = K. We have 
also stressed that the approach unihes static and dynamical properties, generalizing theorems well known for fermion 
systems (ssl - l^ : When the thermodynamic equilibrium is stable, i.e., when the trial free energy f{R} is minimum at 
{R} — the matrix F of second derivatives is positive, and this implies that the eigenvalues of zCF are real. The 

latter property ensures dynamical stability: a small deviation of {R} around is never amplified ('Sec. lIX Cl) . 

The approximate equations (11821) and (I186p determine the expectation values {Qj)^ at the initial time t, and for 

arbitrary t, express these quantities in terms of the expectation values Ra^ ft) = Tr Mq, 2?^°^ ft), which are the symbols 
of the elements Mq, of the Lie algebra. Thus, at first order in the sources, these elements Mq are replaced by the 
c-numbers Ra\t) as if the algebra were replaced by a commutative one. Accordingly, the approximate dynamical 
equations then have a classical Lie-Poisson structure (Sec.|X|. At second order in the sources, we have shown fSec. lVHI) 
that in the initial state the fluctuations and correlations B of the operators Mq have the same form as if the Lie 
algebra were replaced by a mapped algebra in which the commutators [Mq, M^j] are replaced by c-numbers (|109l) . 
which are their expectation values (IM)) . We have also seen that this amounts to replace the deviations p38D by linear 
combinations of operators of non-interacting quasi-bosons (or of harmonic oscillators) together with quasi-scalars. 
This property extends for H = K to two-time correlations. We thus acknowledge two successive modifications of 
the Lie algebra: For hrst-order quantities, the operators Mq are approximated by scalars, whereas for second-order 
quantities it is their commutators which are approximated by scalars. 

Although formal and technical, the present approach is systematic, flexible, and consistent owing to the optimization 
of the generating functional. It allows a variational evaluation of correlation functions, which lie beyond the realm of 
standard mean-held theories. It generalizes some known approaches to arbitrary Lie groups, setting them within a 
unihed framework. It may be applied to various questions of statistical physics and held theory, static or dynamic. 
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at zero or non-zero temperature, at equilibrium or off equilibrium. 
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Appendix A: Derivation of the correlation matrices B and 
1. Ordinary correlations 

As indicated in Sec. lVDl the determination of the correlation matrix IB involves the solution at first order of the 
coupled equations (H1)) . (H5]) for V{t),A{t) in the range 0 < t < /3, with the boundary conditions Da^(O) = 0 and 
A^a\l3) = Mq,. The zeroth order operators (r) = exp[—= exp[(r—/3)K(°)] andD*^°^ = exp[— /3 
where we denote for shorthand have been derived in Sec. lIV Al At first order, we have parametrized 

the combination V{t) A{t) by ~ that is. 


i?«(r) ^ Tr[D(o)(r)Mii)(r) +Dii)(r)M(°)(T)](M;3 - <VTrp(°) . 


The matrix element was then identified with 

M^p = R^% = l3) 


(Al) 

(A2) 


= Tr M„(M, - ^ . 


P ' TrT>(o) 

A preliminary task is to find the r-dependence of (r). This is done through the equation 

^[V{r)A{T)]=[V{T)A{T),K{R^^}] , 

a consequence of (Hill . lHSl) equivalent to 

dRp^jT) 
dr 


= -ihTi^K.'^{R^^} Rf^ir). 


(A3) 


(A4) 


This equation, since the boundary conditions Va'^ (r = 0) = 0 and Aa'^ (t = /3) = are imposed at different times 
T, cannot fully determine the quantity but it will provide its dependence on r. Expanding (IA4I) up to first 


order involves the kernel 


K{R^^{t)} K° + i?W(T) Ms , 


(A5) 


which depends self-consistently on ^(t) and Mq^(t) through {R^^\t)}. Using (IA5|) . then (1551) . (|551) and (|55)) . one 
obtains 


=-ihCi3^ R'd^ (r), 
dr 


which is readily solved, in terms of the still unknown quantities = /3) = B^/s, as 


= 


, i fl C F (/3 —r) 


i?«(/3)=B„^ 


^ 2 /IFC (r — 


(A6) 


(A7) 


We now turn to the determination of Va^ (t) which enters (IA2I1 for t = f3. (We shall not need Aa'^ (r) for 0 < r < /3.) 
The first order contribution to (IT4)) takes the form 


dr 




(A8) 
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Using the boundary condition 'Da\o) = 0, we can solve (IA8|) as 

(t) = - f . 

Jo 


(A9) 


The integrand involves the transform of the element of the Lie algebra by the element ° of the group, 

with . This transform is also an element of the algebra given by the automorphism (1271) for = Pfo) 

and A = (r — r')//3 as 


„K(°>(r'-r) |V/| _ 




n <5 


Mi. 


Inserting (IA7|) and (lAlOl) into (IA9I) leads to 

P«(r)^W(r) = 


(AlO) 

(All) 


- / dr'] 


^a/3 


gZfiFC(T'-/3) _|_ ^-1 gi?!,/3 '-CS(t'- 
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We note that the bracket in Ea. (IAll|) is the derivative with respect to t' of 


gifiFC(r'-/3) _ I 


]Pgi?i/3 ^CS(t'-t) _|_g 


^CS{t'—t) _I 


ih¥C ihCS 

where (e^ — l)/x is defined by continuity as 1 for cc = 0. Thus, by integration of (lAllI) . we obtain for r = /3: 




•-I 


■ —zfl/9FC ][ ^ — ihi 

f-+ - 

ih¥C ihCS 


-I M 




(A12) 


We have thus found the expression of 'Da\p) to be inserted into (|A2p so as to obtain The resulting two 

terms in (IA2I) involve a correlation in the state oc exp[—between two operators of the basis {M} of the 
Lie algebra, Mq, and M/g, or M/g and M.^, respectively. Such correlations have been evaluated in Sec. lIII Cl (within 
replacement of V by and are provided by Eq. dm. Inserting both terms in (IA2I1 . we find 

(A13) 


J _ gifiCS 


-B F 


,-inpcF 


-I 


— i hi 


,-ihSC 


-I 


i?iCF 


insc 


S S' 


-1 


ihsc 

g — i hCS> _J ’ 


and hence 


i?iCF 

_e ?i/3CF 


F- 


(A14) 


2. Kubo correlations 


In order to find a variational approximation for the Kubo correlation 

I _ _ _ 

Tr- / dre'^^ MpD-TiMaDTrUpD, (A15) 

P Jo 

we first replace A{t) by I in the generating functional lnTrA(<i) D and introduce, instead of the sources ij{t) entering 
the exponent of A, small sources 5J°‘ entering terms —added to K in the exponent of I? = exp(—/lAT). We 
now deal with a single exponential operator instead of a product of two, and Kubo correlations are the second-order 
terms in {JJ} in the expansion of InTrZ?. 

The variational approximation B^ for Kubo correlations is then obtained through the formalism of Sec. lIV Al within 
replacement of the image K by K — This yields lnTrexp(—/SAT -I- (5 J“Mq) ~ —/3/{i?} -I- 6J°‘Ra where 

{R} = -I- 6R} is given by the stationarity condition f3df{R^°^ + SR}/dRa = f3¥°‘^6Rp = 6J°‘. As in the 

evaluation of the thermodynamic coefficients, the second-order terms in J} are finally found as 

Ba/3 = ^(F-')a/3 = [(/3K-S)-l]„;3. (A16) 

Here as in the case of B, this approximation for (IA15|) is obtained by replacing —S by /3F = /3K — S in the naive 
expression (1251) written for V = 'D^^'>. 
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Appendix B: The variational results and the effective state 


We prove in this Appendix that the variational approximations Trivial) ~ (Ma)app = Tr MqM^_D ~ 

= Map + Ra^ found in Secs. IIVI and IVl are reproduced in the mapped Hilbert space as exact 

expectation values of Mq and MaM/? over the effective state defined by Eg. (11111) . 

Let us introduce the generating function 


(/>{A} = In 


Tr exp( — /3 F + A“ Ma) 
Tr exp( — /3F) 


(Bl) 


which will produce the expectation values (Mct)map and the Kubo correlations (Ma M/3)jjjap ™ state F by derivation 
with respect to the sources A“. We can rewrite it as 


^{A} = In■ 


Tr exp 


-i/3MaF“T'M;^ + A“i?i°^ + i/3-iA“(F-i)^ A'^ 


Tr exp(-/3£) 

where the operators Ma in the mapped space are defined through the shift (we drop Mo)- 


(B2) 


(B3) 


These operators obey the same commutation relations as (11091) . so that the replacement of {M} by {M } does not 
modify the trace. Hence, we find 


0{A} = A“i?W + ^A“(F-i)^^AW 

We now expand Tr exp[ — /3 F + A“ Ma] in powers of the sources A“: 

"Pp g “ /5 F+-^“ Me ~ 

Tre-'^^ +A“Tr(e-'^^Ma) + ^^“AT' f dr Tr (e" M, 

2 P 7o 

Inserting in m and identifying with (IB4I) . we recover at first order 

(Ma)map = R^°^ ■ 

At second order, we obtain in the space the Kubo correlations of the operators {M}: 

((Ma-^i°^)(M^--Rr))lp = 

1 


(B4) 


(B5) 




(B6) 


(B7) 


i^''dre-^(Ma-i?W)e--^(M^-i?P)\ = ^ (F-')a^ 

/ I 


map 


In order to derive therefrom the ordinary correlations we proceed as in Sec. lHI Cl We note that 

e^-Ma e[F, Ma] ef d(CF)J(M^ - i?®)e“, 


(B8) 


where we used the expression (Hill) of F and the commutation relations H09I) . Integration over r then yields 

e"^(Ma - F®) (M;5 - , (B9) 

and hence, through a new integration as in (l40l) . 


((M., - = 


(BIO) 
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By combining (11121) and (IBIOI) we find that the ordinary correlations of the operators {M} are given by 

((M, - 4»')(M, - 

SO that we recover here the matrix B of Eg. (155)1 . now derived as an exact correlation in the mapped space M’ . 
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